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FOREWORD 


This  document  presents  the  theoretical  results  of  a  study  designed  to  devive  tech¬ 
niques  for  constructing  experimental  mathematical  models  of  linear,  time-invariant, 
network-like  control  systems.  The  study,  which  is  addressed  fundamentally  to  the 
identification  problem  in  control-system  theory,  covers  two  principal  areas  (1)  time- 
domain  techniques  for  determining  system  and  system  input  matrices  and  (2)  the  ex¬ 
traction  of  parameter  submatrices  from  system  input  and  system  matrices.  The  report 
discusses  the  basic  philosophical  approaches  taken  by  others  in  this  field  toward  solving 
the  Identification  problem,  Including  the  model  adaptive  approach,  the  black  box 
approach,  computational  techniques,  quasilineariration,  the  differential  approximation, 
the  algebraic  realitation,  the  instrumental  variable  technique,  and  the  algebraic  mode! 
adaptive  technique.  While  all  these  approaches  are  considered,  the  basic  one  applied 
In  this  study  is  that  of  model  adaptive  Identification. 

This  work,  which  is  directly  applicable  to  various  current  NWC  control-system  pro¬ 
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CHAPTER  I 


INTRODUCTION  AND  BACKGROUND 


INTRODUCTION 

The  general  topic  of  the  thesis  deals  with  what  has 
come  to  be  known  in  the  control  system  literature  as  the 
"identification"  problem.  Given  a  physical  system  which 
is  excited  with  an  input  function  _u  ( t_)  ,  ts(0,T]  ,  the  cor¬ 
responding  output  function  *  te[C,Tj  is  measured.  One 

is  interested  in  using  the  input-output  data  to  construct 
a  mathematical  model  which  will  produce  an  output  which  is 
sufficiently  close  to  the  measured  output  when  the  model 
is  excited  with  the  measured  input. 

There  are,  of  course,  many  variations  of  the  above 
viewpoint.  A  broad  classification  is  into  stochastic  and 
deterministic  identification  problems. 

In  the  stochastic  identification  problem,  one  allows 
that  the  measurements  mvc  noise  superimposed  on  the  true 
input-output  functions.  No  noise  is  assumed  present  m 
the  deterministic  case.  Further,  within  each  of  the  sto¬ 
chastic  and  deterministic  classifications,  there  ■'rc  sub¬ 
class  1  f  i  c  a  1 1  o  n  s  into  linear  and  nonlinear  identi l ication 
problems.  These  subcl assi t ications  refer  to  the  general 
properties  of  the  model.  Each  linear  or  nonlinear 
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identification  problem  also  may  be  categorized  into  model 
types  such  as  static,  differential  equation,  difference 
equation,  transfer  function,  and  pulse  transfer  function 
types  - 

Differential t static  and  difference  equation  models 
may  be  called  either  multi-input  or  multi  output  models  if 
there  are  more  than  one  input  and  output.  Transfer  function 
and  pulse  transfer  function  models  are  generally  single 
input-single  output  models.  Additionally,  differential 
equation  and  transfer  function  types  have  a  continuous  in¬ 
dependent  variable,  generally,  time.  Difference  equation 
and  pulse  transfer  function  types  have  a  discrete  independ¬ 
ent  variable. 

From  the  above  discussion,  there  exists  a  complex 
system  of  classifying  the  types  of  models  which  may  be  used 
to  try  to  match  a  given  set  of  input-output  functions.  No 
guarantee  exists  that  a  given  model  can  match  a  given  input- 
output  function,  hence,  considerable  insight  and  knowledge 
of  the  tested  system  may  be  necessary  to  choose  a  model 
which  will  match  the  measured  input-output  data  within  the 
desired  error. 


BACKGROUND 


A  convenient  method  of  discussing  solutions  to  the 
identification  problem  is  to  classify  some  basic  philosophi¬ 
cal  approaches,  and  then  to  discuss  the  various  computational 
techniques  applicable  to  each  approach. 

The  Model  Adaptive  Approach.  -  In  the  model  adaptive 
approach,  one  assumes  a  model  for  the  system  of  known  struc¬ 
ture.  By  the  term,  known  structure,  it  is  meant  that  every¬ 
thing  is  assumed  known  but  the  parameters  of  the  model. 

That  is,  if  there  are  system  nonlinearities,  the  general 
bri..  of  the  nonlinearity  is  known  but  any  multiplier  or 
additive  values  are  unknown.  Also,  the  form  of  the  matri¬ 
ces  which  describe  the  system  are  known.  A  typical  example 
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would  be  a  companion  matrix  where  at  most  n  of  the  n  pos¬ 
sible  matrix  elements  are  unknown.  After  the  assumption 
of  a  model,  the  identification  problem  then  becomes  one  of 
adjusting  the  model  parameters  in  such  a  fashion  that  the 
error  between  the  model  response  and  system  response  to  the 
same  input  is  sufficiently  small  in  some  sense.  All  avail¬ 
able  knowledge  and  physical  principles  pertinent  to  the 
tested  system  are  useful  in  the  model  adaptive  approach. 

Such  knowledge  includes  constraints  on  the  parameter  values, 
types  of  nonlinearities,  and  system  order. 
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The  Black  Box  Approach.  -  In  the  black  box  approach, 
one  assumes  as  little  as  possible  about  the  tested  system. 
The  usual  expectation  is  that  the  system  is  linear,  but 
test  signals  of  varying  amplitude  are  used  to  confirm  the 
linear  expectation.  Most  commonly,  the  test  signal  used  is 
a  sinusoidal  function  and  the  results  are  typically  dis¬ 
played  as  a  frequency  response  curve  of  amplitude  and  phase 
as  a  function  of  input  frequency.  The  single  input -single 
output  type  of  system  is  the  most  common  target  of  the 
black  box  approach. 

Computational  Techniques.  -  Many  and  varied  types  of 
computational  techniques  have  been  developed  primarily  for 
the  model  adaptive  approach.  A  brief  discussion  of  the 
dominant  techniques  is  presented. 

Quasi linearization.  -  Quasilinearization  is  a  compu¬ 
tational  technique  developed  primarily  by  Bellman  and  his 
coworkers  [IT-1] .  Extensions  of  the  technique  have  been 
made  by  Lee  [IT  4]  and  Sage  [CT-4] .  It  is  an  iterative 
procedure  based  on  a  generalization  of  the  Newton -Raphson 
technique  for  finding  the  roots  of  a  nonlinear  equation. 
Both  linear  and  nonlinear  systems  are  modelled  using  quasi - 
linearization.  If  the  technique  converges,  the  convergence 
is  quadratic  in  nature.  There  is,  of  course,  no  a  priori 
guarantee  of  convergence.  Also,  the  region  of  convergence 
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may  be  large  or  small.  Measurements  required  are  the  system 
input-output  measurements. 

Differential  Approximation.  -  Differential  approxima¬ 
tion  is  a  computational  technique  where  knowledge  of  the 
system  inputs,  states  and  state  derivates  is  required. 
Bellman  [IT-1]  and  his  coworkers  were  instrumental  in  de¬ 
veloping  this  technique.  It  is  basically  a  gradient  tech¬ 
nique  where  an  integral  of  a  suitable  norm  of  the  differ  - 
,‘nce  between  the  known  state  derivative  and  the  analytical 
expression  (dependent  upon  known  states  and  unknown  parame 
ters)  for  tin  state  derivative  is  minimized.  Minimization 
of  the  integral  yields  a  set  of  nonlinear  algebraic  equa¬ 
tions  in  the  parameters.  It  is  often  used  as  a  starting 
routine  for  other  more  sophisticated  techniques.  Linear 
and  nonlinear  systems  may  be  modeled  by  this  technique. 

Algebraic  Realization.  -  Identification  using  alge 
braic  realization  theory  was  developed  by  one  of  Kalman's 
students  at  Stanford  University  [IT-3] .  The  technique  is 
based  on  the  concept  of  detcmining  the  Markov  parameters 
of  the  system.  The  Markov  parameters  are  an  infinite  se¬ 
quence  of  matrices  which  are  related  to  the  system  parame¬ 
ters.  It  is  assumed  that  the  system  is  linear  and  time 
invariant.  Based  on  the  above  assumption,  the  identifica¬ 
tion  problem  is  composed  cf  two  sub-problems.  The  first 


sub-problem  is  concerned  with  linear  algebraic  operations 
on  the  data  to  determine  a  finite  number  of  the  Markov 
parameters.  Secondly,  linear  algebraic  operations  on  the 
Markov  parameters  are  used  to  determine  output,  system  and 
input  matrices  of  the  system.  There  is  a  certain  amount 
of  arbitrariness  in  the  second  part  of  the  problem.  Conse  • 
quently,  the  results  are  not  unique,  but,  vary  within  a 
similarity  transformation.  The  required  measurements  are 
input  and  output  data. 

Instrumental  Variable  Technique.  -  Young  [IT-10, 

IT-11]  and  Wong  [IT-12]  have  developed  computational  pro¬ 
cedures  using  instrumental  variables.  The  technique  is  de¬ 
signed  primarily  for  single  input -single  output  linear, 
time -invariant  systems.  Noisy  input -output  data  is  meas¬ 
ured  and  used  in  recursive  algorithms  for  real  time  parame¬ 
ter  estimation.  To  bypass  the  necessity  to  determine  de¬ 
rivatives  of  the  input  and  output  explicitly,  the  input 
and  output  are  passed  through  a  series  of  "state  variable" 
filters  whose  outputs  constitute  the  instrumental  variables. 
By  proper  choice  of  the  filters,  the  instrumental  variables 
and  parameters  have  the  same  relationship  to  each  other  as 
do  the  original  system  states  and  parameters.  Consequently, 
by  measuring  the  instrumental  variables,  the  parameters  can 
be  computed. 


6 


Algebraic  Model  Adaptive  Technique.  -  Yore [IT-9]  pre¬ 
sented  a  thesis  which  considered  two  separate  but  related 
problems.  The  first  problem  was  determining  the  system  and 
input  matrices  of  a  linear  time-invariant  system  where  the 
entire  state  and  input  vectors  were  measured.  An  iterative 
procedure  based  on  linear  algebraic  techniques  was  used  to 
cause  a  model  state  to  converge  to  the  measured  state.  The 
second  problem  was  to  assume  that  the  system  generating  the 
data  was  a  network-like  system,  and  to  devise  methods  of 
determining  the  values  of  the  components  contained  in  the 
system  and  input  matrices.  For  example,  if  a^  =  rlcl 
an  element  of  the  system  matrix,  what  is  the  value  of  r^c^? 
Yore's  work  was  apparently  the  first  to  address  the  so- 
called  "component  parameter"  identification  for  network - 
like  systems.  The  "component  parameter"  identification  is 
the  second  part  of  the  present  thesis. 
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PROBLEM  STATEMENT  AND  ASSUMPTIONS 


Yore's  [IT-9]  work  formed  the  motivation  for  the 
present  thesis.  It  was  felt  that  his  techniques  were  use 
ful  but  could  be  improved. 

Basic  Assumptions.  -  The  following  assumptions  were 

made: 

1.  The  measured  data  is  noiseless  and  generated  by  a  net- 
work-like  system. 

2.  System  topology  is  known,  i_ .£. ,  it  is  known  what  types 
of  components  are  present  and  how  they  are  interconnected. 

3.  Actual  parameter  values  are  unknown. 

4.  The  system  is  linear  and  time  invariant  over  the  meas¬ 
urement  interval. 

Problem  Statement.  -  Determine  techniques  and  condi¬ 
tions  to  identify  the  parameters  of  the  system  described 
in  the  above  assumptions. 
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CHAPTER  II 

RESULTS,  CONCLUSIONS,  AND  RECOMMENDATIONS 


INTRODUCTION 

A  summary  of  results  and  comparisons  with  some  cor¬ 
responding  results  of  Yore  is  presented.  Also  included  in 
the  chapter  are  conclusions  based  on  the  results  and  recom¬ 
mendations  for  further  research. 


SUMMARY  OF  RESULTS 

A  chapter-by-chapter  summary  of  the  research  results 
is  given. 

Chapter  III,  Model  Adaptive  Identification.  -  The 
algorithm  which  adapts  the  model  to  fit  the  measured  state 
is  based  on  linear  algebraic  computations.  Both  input  and 
system  matrices  can  be  identified  using  measurements  of  the 
state  and  normal  operating  inputs.  The  algorithm  is  itera¬ 
tive  and  proceeds  from  iteration- to-iteration  by  the  fol¬ 
lowing  relations. 


(1) 


.  n  +  1 
vm 


m 


+  A  A 


n 


(2)  an+1  =  ♦  a  b 

m  m 


n 
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The  matrices  are  the  model  system  and  input  matrices 


at  the  n-th  iteration,  respectively.  Computation  of  the 
updating  matrices  is  accomplished  by  the  following  discrete 
time  equivalent  of  the  continuous  time  error  differential 
equation. 


(3)  _an  (k+lT)  =  e 


AmT  f  n 


{en(Kt,+(  )} 


The  error,  e.n™  x  -  x^#  is  defined  by  the  measured  state,  x# 


and  the  computed  model  state  x”.  In  (3)  u  is  the  input 


and  T  is  the  sampling  interval.  By  taking  a  sufficient 
number  of  samples,  AAn,  ABn  are  computed  using  matrices 
whose  columns  are  defined  by  (3) . 

A  problem  which  could  occur  is  that  the  model  has  been 
chosen  of  incorrect  order.  Two  possible  effects  of  incor¬ 
rect  order  are  that,  (1)  the  algorithm  does  not  converge, 
and  (2)  the  algorithm  converges  to  erroneous  A, B.  The 
important  consideration  here  is  which  of  the  states  are 
neglected.  By  neglecting  strongly  excited  states,  poor 
results  can  be  achieved. 

Another  problem  which  could  occur  is  that  the  state 
and  input  could  be  measured  with  some  error.  If  measure¬ 


ment  error  occurs,  then  the  algorithm  will  converge  to 
erroneous  A,B  matrices,  if  it  converges  at  all.  The  magni¬ 
tudes  of  these  errors  depend  on  the  A,B,u,  and  of  course. 


WWJPWWWW^ 


1  isr^JwwwiT^i  w 


'? 

i 

the  measurement  error.  The  sample  interval  and  norm  of  the  j 

i 

system  matrix  are  related  to  the  error  in  identification  at  } 

the  first  iteration.  Nearly  one  step  convergence  can  be  ' 

achieved  under  certain  conditions.  Some  numerical  experi¬ 
ments  for  the  first  order  free  case  indicate  that  converg¬ 
ence  is  guaranteed  for  a°  -  0  if  the  following  conditions 
exist: 

(4)  | aT  j  -  arbitrary  a<:  0 

(5)  aT  <  2.5  a  >  0 

Example  3  showed  that  under  certain  conditions,  the 
requirement  to  measure  the  entire  state  vector  could  be 
relaxed  by  identifying  a  model  whose  state  could  be  re 
la ted  to  the  measured  output. 

Chapter  IV,  Parameter  Estimation  Using  Nonlinear 
Observers.  -  Adopting  the  viewpoint  that  the  system  parame¬ 
ter  could  be  collected  into  a  vector,  £,  with  differential 
equation,  £  =  o,  the  original  linear  system  could  be  aug¬ 
mented  to  form  a  nonlinear  system.  Using  the  Kronecker 
product,  (6)  through  (9)  show  the  form  of  the  augmented 
system. 

(6;  x  =  Ifiu^a  +  I&uSd 

(7)  a  =  o 

(8)  b  =  o 
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(9)  2  =  Cx 


The  vectors  jt,  b  are  formed  by  collecting  the  rows  of 
the  system  matrix  A  and  the  input  matrix  B  into  the  vectors 
_a,b,  respectively.  Not  all  of  the  states  are  measured.  The 
output  matrix  C  is  of  order  mxn,  and  the  input,  system 

matrices  are  of  order  nxp,  nxn,  respectively. 

2 

A  nonlinear  observer  of  dimension  n+n  +np  is  con¬ 
structed  whose  state  vector  asymptotically  approaches  the 
vector  (x,a,b)t  under  certain  conditions.  Success  of  the 
scheme  depends  on  one's  ability  to  stabilize  the  observer. 
Since  the  observer  is  nonlinear  and,  generally,  time-varying, 
stabilization  may  be  difficult.  A  procedure  for  stabiliza¬ 
tion  was  developed. 

Chapter  V  -  Estimation  of  Parameter  Sub-Matrices. 
Chapter  v  assumed  that  the  system  to  be  identified  was  a 
network-likc  system.  Therefore,  the  A,B  matrices  have  a 
special  structure.  Assumption  of  the  network-like  system 
constrains  the  elements  of  A, 8.  The  constraints  appear  by 
forcing  certain  elements  of  A,B  to  be  either  zero  or  one. 
Other  elements  are  related  by  algebraic  equations  to  each 
other.  Thus,  information  about  A,B  can  be  determined  prior 
to  actual  testing  of  the  system. 
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Knowledge  of  the  state,  input  vectors,  and  A,B  matri¬ 
ces  may  not  be  sufficient  to  determine  the  parameter  matri¬ 
ces  uniquely.  Some  results  are  given  which  show  where  ad¬ 
ditional  measurements  should  be  taken  to  give  unique  results. 
Both  linear  and  nonlinear  techniques  for  determining  the 
system  parameters  are  given.  It  is  shown  that  for  a  special 
class  of  networks,  all  parameters  may  be  determined  from 
the  A, B  matrices. 

COMPARISON  WITH  YORK'S  RESULTS 

The  notation  in  Chapter  ffi  is  used  to  compare  the  re¬ 
sults  achieved  here  with  those  by  Yore.  Table  1  shows  the 
difference  in  the  algorithms. 

It  is  the  author ‘s  opinion  that  Yore  unnecessarily 

restricted  his  technique  by  requiring  the  sample  interval 

AT 

be  sufficiently  small  that  e  •»  i+a?. 

If  Yore's  time  requirement  is  used  for  the  author's 
algorithm,  then  one  step  convergence  is  achieved.  One  step 
convergence  is  not  achieved  for  the  same  sample  time  using 
Yore's  algorithm. 

Two  restrictions  on  the  sample  time  are  given  by  Yore. 
For  the  first  order  autonomous  case,  it  is  necessary  and 
sufficient  that  jaTj  <_  Ln2  for  convergence  of  the  algorithm. 
For  higher  order,  non  diagonal,  autonomous  case,  it  is 
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sufficient  that  n  max 

i#  3 


Ln2  for  the  identification 


to  converge.  The  results  of  Chapter  III  show  that,  at 
least  for  the  first  order  case,  the  sampling  time  require¬ 
ments  are  much  less  restrictive  than  Yore  presents. 

Comparing  Yore's  results  with  those  in  Chapter  V  on 
determination  of  the  parameters  in  the  a,B  matrices  is  dif 
ficult  because  the  approaches  are  fundamentally  different. 

Yore's  main  effort  was  characterization  of  a  suffi¬ 
cient  number  of  measurements  in  terms  of  his  concepts  of 
statically  and  dynamically  independent  variables.  His  re¬ 
sults  are  phrased  in  terms  of  the  existence  of.  such  sets  o 
variables.  No  procedures  or  techniques  are  given  to  de¬ 
termine  which  network  currents  and  voltages  constitute 
statically  and  dynamically  independent  sets  of  variables. 
No  direct  use  of  system  topology  was  made. 

In  contrast,  the  effort  in  Chapter  V  was  based  ex¬ 
clusively  on  detailed  knowledge  of  system  topology.  The 
main  effort  was  to  use  the  topology,  A,B  matrices,  state, 
and  input  vectors  to  determine  the  parameter  submatrices. 

As  a  byproduct  of  the  computations,  the  need  for 
additional  measurements  other  than  state  and  input  rs  de¬ 
termined.  Direct  application  of  linear  algebraic  tech¬ 
niques,  principally,  determination  of  ranks  of  various 
matrices,  is  required. 
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CONCLUSIONS 


The  basic  algorithm  in  Chapter  III  is  effective  and 
works  well  for  the  cases  in  which  it  is  applicable.  Its 
capability  of  providing  nearly  one  step  convergence  should 
prove  useful.  The  two  principle  effects  of  measurement 
no lie  anu/or  incorrect  model  order  are,  (1)  non-uniqueness 
of  results,  and  (2)  unique  but  erroneous  results. 

Measurement  noise  can  be  significant  if  the  states 
are  weakly  excited.  The  data  matrix  which  must  be  in¬ 
verted  becomes  ill-conditioned.  A  consequence  of  the  ill- 
con  iitioning  is  that  small  errors  in  the  data  matrix  ele¬ 
ments  cause  large  errors  in  the  updating  matrices.  The 
errors  can  cause  the  iterations  to  oscillate  about  the 
correct  values  of  A, B  matrices.  A  filter  can  be  cascaded 
with  the  basic  algorithm  to  smooth  the  fluctuations. 

Knowledge  of  the  system  structure  can  be  used  to  de¬ 
velop  a  nonlinear  observer  which  has  on-line  potential  for 
combined  state  and  parameter  estimation.  The  observer  is 
best  suited  for  analog  or  hybrid  mechanization  because  a 
large  number  of  first  order  differential  equations  with 
relatively  widely  separated  time  constants  need  to  be 
integrated. 

The  critical  problem  is  insuring  that  the  observer  is 
stable.  A  complicating  factor  is  the  requirement  for  the 
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feedback  for  each  parameter  equation  to  decrease  to  zero  f 

\ 

while  keeping  the  overall  system  stable.  A  technique  was  j 

i 

developed  to  stabilize  the  observer.  j 

Knowledge  z£  the  system  topology  is  necessary  for  s 

unique  parameter  determination  for  the  network-like  case.  j 

Straightforward  results  and  techniques  have  been  demon¬ 
strated  for  parameter  determination.  | 
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. . .  ......... 


RECOMMENDATIONS  FOR  FUTURE  WORK 


Future  research  effort  should  be  directed  to  answer 

the  following  questions.  Most  work  should  be  directed  to 

the  computational  aspects  of  the  problem. 

1.  How  can  the  effects  of  widely  spread  eigenvalues  be 
handled  computationally? 

2.  If  the  measurements  are  corrupted  with  noise,  how 
can  the  results  of  statistical  estimation  theory  be 
applied  to  give  "good"  estimates  of  A, B  using  the 
algorithm  of  Chapter  3? 

3.  How  can  the  effects  of  incorrect  model  order  be  mini¬ 
mized? 

4.  How  can  the  sample  interval  for  the  basic  algorithm 
be  chosen  for  most  rapid  convergence?  Should  it  be 
modified  from  iteration  to  iteration? 

5.  How  can  the  model  adaptive  algorithm  be  used  where 
only  a  subvector  of  the  state  vector  is  measured? 

That  is,  when  can  additional  operations  on  output 
data  such  as  integration  be  used  to  identify  another 
model  whose  parameters  are  related  to  the  original 
parameters  by  known  relations? 
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CHAPTER  HI 


MODEL  ADAPTIVE  IDENTIFICATION 


INTRODUCTION 

A  model  adaptive  algorithm  is  derived  under  the  usual 
assumption  that  the  pair  [A,B]  is  completely  controllable 
[Kalman,  CT-2] .  Further,  it  is  assumed  that  the  state 
vector  and  input  vector  are  measured  perfectly.  After  the 
algorithm  is  developed,  the  effects  of  measurement  error 
and  incorrect  model  order  are  considered.  Some  uniqueness 
results  and  computational  rules  of  thumb  are  presented. 

The  physical  system  is  assumed  to  be  governed  by  the 
constant  coefficient,  vector-matrix  differential  equation 
(10). 

(10)  _x(t)  =  Ax  +  Bvi(t)  ,  x^  =  x. (0)  te[0,tf] 

and  that  x(t)  and  ujt)  are  n  and  p  vectors  measured  per¬ 
fectly  on  the  interval  [0,t^].  A,  B  are  not  known. 

A  model  of  the  system  is  assumed  to  have  the  constant 
coefficient  vector-matrix  differential  equation  (11) . 

UD  4  =  V,  +  V'  xm(0)£  4 

The  chart  records  of  x(t)  *  u.(t)  are  divided  into  n+p+1 
samples  which  are  used  to  determine  the  A,  B  matrices. 
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as  . - 


By  assuming  a  pair  £Am*Bm],  integrating  (11)  and  com¬ 
paring  the  measured  and  model  states  at  the  sampling  instants 
t^,i  =  0, 1, 2, . . . , n+p+1,  the  model  is  forced  to  converge  to 
the  system  represented  by  (10) .  The  updating  scheme  is 
algebraic. 

DERIVATION  OF  ALGORITHM 

The  method  of  forcing  the  model  data  to  converge  to 
the  measured  data  is  to  define  the  following  quantities. 

AA,  AB  are  the  updating  matrices. 

(12)  A  =  Am  +  AA 

A 

(13)  B  =  Bm  +  AB 

In  terms  of  the  defined  variables  in  (12)  and  (13) , 
equation  (10)  is  rewritten. 

(14)  x  -  A  x  +  AAx  +  B  u  +  ABu 

—  m  m~“  — 

Let  e(t)  =  x  -x  (t)  denote  the  difference  between  the 

~  —  “TO 

measured  data  and  the  model  data,  then  by  subtracting  (11) 

from  (14)  the  error  differential  equation  is  (15) . 

A 

(15)  e  =  A  e  +  AAx  +  ABu,  e(0)  =0 

—  ur-  *—  —  —  — 

Since  e,  x,  u,  are  known,  AA,  AB  can  be  computed.  In 
order  to  compute  AA,  AB  (15)  is  converted  to  the  equivalent 
discrete  time  system  (16) . 
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(16)  e(k+l) 


=  eAmT  {  e(k)  +  (QJ  e  MO  (AAAB)  (u(k)  )} 


where  T  =  t^+^-t^,  i  =  0, 1.2, . . .,n+p  and  is  the  sampling 
interval. 

The  error  and  input-output  vectors  are  formed  into 
the  following  matrices. 

UO  =  ^u(o),...,  u(n+p-lj) 

XO  =  (x (o) , . . . ,  x(n+p-lj) 

EO  =  (e  (o) , . . . ,  e  (n+p-1)) 

El  =  (.e(l) , e^(n+p)  ) 

Using  (16),  the  matrices  defined  by  (17)  -  (20)  are 
related  as  follows: 

(21)  El  =  e  m  {  EO  +  (QJ  e  Mr)  (AaAB)  (uo)} 

Assuming  the  existence  of  the  inverse  of  (^q),  then» 

the  updating  equations  to  improve  the  model  in  an  iterative 
fashion  are: 


(17) 

(18) 

(19) 

(20) 


<22>  \k+1  -  C +  *a 


(23)  Bmk+1  =  B  k  +  LB 
m  m 
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(24) 


(AAAB)  = 


(JWX  {.^V* 


-l 


EFFECTS  OF  MEASUREMENT  ERROR 

The  effects  of  measurement  error  on  the  A,B  matrix 
errors  are  formulated  in  the  present  section. 

Assume  that  the  actual  state  and  input  vectors  are 
measured  with  errors  e  and  e  and  define  the  errors  by 

* — “X  •—‘U 

(25)  -  (28). 

A 

(25)  x  =  x  +  e 

"  "“A 

(25)  u  =  u  +  e 
—  —  — u 

A 

(27)  A  =  A  +  EA 
A 

(28)  B  =  B  +  EB 

EA,  EB  are  identification  errors  induced  by  e  ,  e  . 
The  identification  algorithm  operates  on  the  apparent 

system  given  by  (29) . 

A  A  A  a  A  A  A 

(29)  x  =  A  x  +  B  u,  x(0)  =  x 

Subtracting  the  actual  system  equation  x  =  Ax  +  Bu_,  2i(°)  = 
x^  from  (29)  yields  the  equation  relating  the  measurement 

error  to  the  induced  identification  errors  EA,  EB  in  (30) . 

.  A  A 

(30)  e  =  Ae  +  Be,  +  (EA)x  +  (EB)u 
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Estimating  EA,  EB  using  the  discrete  time  equivalent 


i  ' 


5 


of  (30)  yields  (31)  which  is  the  relationship  of  the 
measurement  error  and  the  A,B  error. 

,T 

— u 


(31)  ex(k+l)  =  e 


AT  J^OO  +  (QI  e  dT)(fi^(k)  + 


(EA  EB) 


x(k)  \ 


u(k) 


thus,  the  errors  are  given  by  (32) . 


L ' 


/  »T  _  \  —  1  i—l 

oJ  e  <*<]  le-^-EQ^  -BBoJj  A,j 


(32)  (EA  EB)  = 

El„.  EO  ,  £q.  &)  are  defined  as  in  the  section  on 


x  x 
the  basic  algorithm. 


EFFECT  OF  SYSTEM  ORDER  ERRORS 

The  present  section  provides  the  formulation  of  the 
effect  of  errors  in  selecting  the  model  system  order. 

Assume  the  system  is  precisely  modeled  by  the  state 
equation. 


1 

[41 

An  A12^ 

fii  ^  | 

(33) 

1 

[-2 

=3  1 

\A21  A22 

—2  i  ( 

>1  \ 
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2  ; 


where  x^  is  a  r-vector  of  the  assumed  model  order  and  is 


an  n-r  vector  of  the  ignored  variables.  Now  letting  e^ 
52  *.]_  and,  B^  =  Bra  +  AB,  A^  =  as 
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basic  algorithm,  (33)  is  transformed  into  a  vector  matrix 
differential  equation  in  the  error  variable  and  the  sub¬ 
vector  of  the  ignored  variables  and  is  given  in  (34) . 


(34) 


<2 

\  Z-2 


/  A 

1  _  m 

A  ' 
A12; 
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A22  j 
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jA21  I 
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The  effect  of  the  ignored  variables  is  that  of  an 
additional  input  and  can  be  considered  as  an  unknown  dis¬ 
turbance. 

Analysis  of  (34)  indicates  that  the  eigenvalues  of 
the  error  plus  ignored  variable  system  matrix  are  those 
of  the  model  system  matrix  and  those  of  the  ignored 
system  matrix  A^* 


An  estimate  of  the  effect  of  the  error  in  system 
order  on  the  identification  can  be  obtained  by  the  discrete 
time  equivalent  of  (34)  . 

Decomposing  (34)  into  two  equations. 


A22T 


(35)  x2(k+l)=e  **  <| 


f  .T 

*2  00 jo'  < 


-A 


\t 

22  r 

d  t 


A21— 1  (k)+B2-(k^  I 


A^T  | 

(36)  £(k+l)«e  <c(k) 

r 


i  _Atn 

+  QJ  o  iAnil  (kjtiBjUlk) 


+h12x2 (k) 
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and  writing  the  matrix  equation. 


T 


(37)  El  =  e 


A  T 
rn 


EO  + 

\ 


(.1  .-\0 


dr' 


'X03 

uo. 


and  solving  for  (AA-^AB^) 
the  case  where  El  =  EO  = 


,T 


+  .V  (J  e-\V  ,x2, 


,  the  induced  error  matrices. 


(0) ,  i.e. ,  x  -»  x., 
— m  —1 


for 


(38)  (ta^)  =  [(QJ  e  l(-fl12X2)]  Q)-1 

and  the  matrix  X2  is  the  n-r  xn+p  matrix  formed  from  n+p 
samples  of  the  solution  of  (35) . 

By  examining  (38),  one  sees  that  if  the  data  matrix 

I  xoi  )  . 

'  uo  ls  non- singular,  then  the  identification  converges 

to  some  (unique)  matrices  A^,  B^  which,  in  general,  are  not 

the  same  as  the  submatrices  A11#  B^.  If  the  data  matrix 

is  singular,  then,  the  identification  diverges  or  at  best 

yields  some  arbitrary  Am»  Bm  depending  on  whether  or  not 

one  desires  to  specify  the  dependent  elements  of  A  ,  B  . 

m  m 

The  important  result  is  that  incorrect  model  order  is 
not  too  critical  to  convergence  of  the  identification. 
Additionally,  identification  using  different  inputs  can 
cause  different  error  matrices  AA^^AB.  The  difference  is 
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caused  by  the  interaction  of  measured  states,  ignored 
states  and  input.  For  example,  if  the  input  did  not  excite 
the  ignored  variables  strongly,  then,  x2(t)  «  o  and  dif¬ 
ferent  (nearly  null)  AA^,  AB,  matrices  will  result  than  if 
the  ignored  states  are  strongly  excited.  A  method  of  check¬ 
ing  for  correct  system  order  of  the  model  could  be  exciting 
the  system  with  different  inputs  and  comparison  of  the  re¬ 
sults.  No  change  in  A^,  with  different  inputs  would 

mean  correct  order.  A  change  in  A  ,  B  with  different  in- 

m  m 

puts  would  mean  incorrect  order.  A  change  in  Am»  with 
different  inputs  could  mean  wrong  order,  but  it  could  also 
be  the  result  of  another  problem. 

It  is  possible  that  the  model  could  be  the  correct 
order  but  some  weakly  excited  states  and  numerical  diffi¬ 
culties  cause  some  changes  in  A  ,  B  using  different  in- 

m  m 

puts.  The  problem  lies  primarily  with  the  computational 
aspects  of  the  identification.  Ill-conditioned  data 
matrices  caused  by  weakly  excited  states  are  difficult  to 
handle.  None  of  the  example  problems  worked  presented  any 
difficulties  of  this  nature. 
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INVERSION  OF  DATA  MATRICES 


A  series  of  results  concerning  the  conditions  for  the 
inversion  of  the  data  matrices  in  the  basic  algorithm  are 
given  under  the  assumption  that  the  data  are  measured  with¬ 
out  error. 

The  following  notation  is  used: 


where  J  is  an  n+px  n+p  matrix  and  X^,  U2  are  nxn  and 
pxp  matrices. 

Theorem  #1 .  Necessary  and  sufficient  conditions  for 
the  inverse  of  (39)  to  exist  are, 

(1)  det  (XitX1  +  ujt^)  ^  0 

(2)  det  (X^X2  +  U^U2)  ji  0 

Proof :  Since  the  inverse  of  ^  J  if  and  only 

if  det  ^  0,  if  a  bound  on  the  determinant  of  (39)  can 

be  established,  then,  the  theorem  is  proved.  Let  A  - 

(A^A2)  33  1  y2j  ,  then,  A^A^  =  and  A^A2  » 

1  2 

*2*2  +  U2U2-  By  Wegner's  theorem  [Bodewig,  M-l,  p.  71), 

dot  h  =  (dot  A^A^}  *  (det  A^A2)  •  Therefore,  det  A  ^  0  if 
and  only  if  the  theorem  is  true. 


Theorem  #2.  Sufficient  conditions  for  (39)  to  have 


an  inverse  are, 

(1)  det  X1  ?  0 

(2)  det  U2  /  0 

(3)  det  (X1-X2U2“1U1)  ^  0 

(4)  det  (U2-U1X1"1X2)  j i  0 


CONVERGENCE  AND  UNIQUENESS 


Theorem  #3.  Sufficient  conditions  for  the  unique¬ 
ness  of  the  identification  are: 

(1)  e_(t)  A  0_  t  e  [0,tf] 

(2)  that  there  exist  n+p  linearly  independent 
vectors  (^[^| )  in  the  open  interval  te(0,tf). 

Proof:  The  error  equation  (15)  is  £(t)  =  A^eft) 

+  (A-Am)j<(t)  +  (B-Bm)u_(t),  _e(0)  £  ()  and  if  _e(t)  A  0  on  the 
closed  interval  te[0,tf],  then  _e(t)  =  £  in  the  open  inter¬ 
val  te(0#tf).  Consequently,  the  error  equation  takes  the 
form  0_  *>  (A-A^I B-B^)  ♦  \~^  )  te(0,tf)  under  the  hypothesis 
of  the  theorem. 

Let  Y  be  the  matrix  of  the  linearly  independent 
samples,  then,  the  matrix  equation  below  is  valid; 

(0)  =  (A-Am:B-B.n)Y 

Since  Y  is  non-singular,  then,  its  multiplier  must 
be  tl;e  null  matrix,  i .  e. ,  A  =  A^  and  B  =  D^. 

Some  insight  into  the  selection  of  the  time  inter¬ 
val  for  the  basic  identification  algorithm  can  be  obtained 
for  the  free  response  case.  The  more  general  case  is  con¬ 
siderably  more  complex  and  has  not  yielded  to  analysis. 
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Theorem  #4.  Let  the  updating  algorithm  (21)  have 

the  initial  value  A  =  (0)  and  consequently  the  first  itera- 

o 

tion,  be  given  by  A^  =  |  (El-EO)XO  then, 

Lim  (A. -A)  =  (0) 

T  -»  0  1 

Proof:  The  proof  proceeds  by  direct  computation. 
Substituting  for  El,  E0,  since  A  ^(0),  the  expression  for 
A1  is  , 

(40)  Ax  =  ~  (Xl-X0m-(XO-X0m)  )  X0"1 


As  a  remark,  (41)  shows  pxecisely  why  one  step  con¬ 
vergence  is  not  obtained,  but  also  shows  that  for  T  suf¬ 
ficiently  small,  that  A^  can  be  very  close  to  A.  A  bound 
on  the  difference  A^-A  can  be  obtained,  however,  the  bound 


will  be  conservative. 


1  AT 

Theorem  #7.  Let  A-^  =  ^  (e  -I) ,  then,  the  error  A-^-A 


is  given  by, 

(43)  ||A,  -A  ||  <  i|A|i2T  e^A^T  =  e,  T,  s>0 

1  ^ 


AT 

Proof:  Using  the  series  expansion  for  e'  ,  it  fol¬ 
lows  that  [Pipes,  M-7,  p.133], 

.  .  1  (a2t2  +  a3t3  +  ...)  a2t  at  ©n  .  n  r  . 

(44)  A^-A  -  7p  .  ~2 —  ~~r~  =  — •  e  *  0  <  0  £  1 


To  get  an  upper  bound,  let  0=1,  and  take  the  norm 
of  (44) . 

(45)  Hawaii  <  lle^ll  <  M2t  JlAllT  =  e 

Equation  (45)  can  be  used  as  a  guide  for  choosing  T 
to  get  within  a  specified  error  on  the  first  iteration. 

Let  €/J|a!I  be  the  desired  fraction  error,  then,  solving  (46) 
for  the  number  ||a||t  and  estimating  ||a||  to  get  an  estimate 
of  T  which  will  achieve  the  desired  maximum  error  fraction 

c/UaU. 

(46)  ||a{|t  ellAHT  =  2  e/  11a  !i 

Figure  (  /  )  shows  [|a||T  vs  e/!|A|| 
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FIGURE 


A  proof  that  a  specific  value  of  aT  which  guarantees 
convergence  for  the  first  order  case  (the  simple  one)  has 
not  been  achieved.  Figures  2  to  4  show  the  results  of 
some  numerical  experiments  for  various  values  of  aT.  Fig¬ 
ure  2  indicates  monotone  convergence  if  aT<0  and  the  possi- 

/ 

bility  of  sustained  oscillation  if  aT>2.5.  Figures  3  and  4 
verify  that  the  convergence  is  monotone  if  aT<0  and  oscil¬ 
latory  convergence  for  aTX).  Based  on  the  numerical  re¬ 
sults,  one  should  use  the  following  guide  for  selection  of 
the  sampling  interval. 

Rule  1.  a<0 

T  can  be  arbitrary  but  should  be  small  to  yield  to 
rapid  convergence. 

Rule  2.  a>0 

T<2.5/a  but  should  also  be  small  to  yield  rapid  con¬ 
vergence. 


The  following  expression  for  the  errors  can  be  de¬ 
rived  by  simple  manipulation  of  the  basic  algorithm. 


(47) 


a  , .  -a 
n+1  n 


n 


aT  anT 
e  -e 

a  T  . 
e  n  -1 


,  n  >  0 


(48) 


i  aT  anT  (eaT-1) 

n+1  ,  aT  t  x  e  -e  n 

=  (e  -1)  -7TJ — ^ -  ,  n  >  0 


n 


a  T.  aT  .. 
e  n  (e  ~X)-1 
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Itera 


CONVERGENCE  PROPERTIES  FOR 
FIRST  ORDER  SYSTEM 

x  =  ax,  a  >  0 
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NUMERICAL  EXAMPLES 


Two  numerical  examples  are  presented  to  illustrate 
the  technique.  In  each  case,  the  initial  assumption  was 
that  A^  =  (0) ,  Bm  =  (0) .  Both  single  input  and  multi-input 
systems  were  identified. 

Example  #1  -  Second  order  1  input  case:  The  first 
example  problem  attempted  was  an  oscillatory  system  de¬ 
scribed  by  the  phase  variable  model. 

<49>  A  =  (-4  -2) 

(50)  b  =  (J) 

A  step  input  of  2  was  used  with  initial  condition 
*^(0)  =  .278,  ^(O)  =  *536.  Runs  were  made  first  assuming 
b  was  known  and  then  assuming  b  unknown.  Convergence  to 
the  correct  A,]3  was  rapid  in  both  cases.  A  sampling  time 
of  .1  second  was  used. 


Table  2.  Identi fication  of  Second  Order  Phase  Variable  System 
With  Known  Input  Matrix. 
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Example  #2  -  Sixth  order,  2  input  complete  circuit:  A 


complete  circuit  in  the  sense  of  Brayton  and  Moser  [PE-3] 
was  chosen  for  the  second  example.  The  network  schematic  is 
shown  in  Figure  5.  Inputs  to  the  system  were  u^  =  20.0, 
u2  =  10.0  +  5  sin  .08  t  and  the  initial  conditions  were  zero. 
The  sampling  interval  was  0.1  sec. 

Again,  it  was  assumed  that  the  A,B  matrices  were  null 
initially.  The  numerical  results  show  rapid  convergence  to 
the  correct  values.  A  general  structure  for  the  A, B  matri¬ 
ces  is  shown  below. 


\ 


f 

! 

1 

r  -  1 

Rici 

0 

1/C1 

0 

0 

(51) 

0  R5=2 

1 

0 

0 

1/C  2 

0 

0  0 

0 

-l/c3 

"1/C  3 

-1/C3 

A1" 

.  -1/LX  0 

1/L1 

"Rl/Ll 

0 

0 

o  -i/l2 

t 

-i/l2 

0 

-r2/l2 

0 

1  0  1/L3 

1/L3 

0 

0 

"VS 

(52) 


B= 


/  0 
:  0 

i  0 


0 

0 

0 

1/L 

0 
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Table  4.  Identification  of  Sixth  Order,  Two  Input 
_ Complete  Circuit _ 


PARAMETER 

0 

1 

4 

ITERATION 

10 

NUMBER 

CORRECT 

VALUE 

A11 

0.0 

-6.15 

-5.27 

-5.01 

-5.0 

A12 

0.0 

-2.30 

.92 

+  .03 

0,0 

A13 

0.0 

10.49 

-4.11 

-  .11 

0.0 

A14 

0.0 

5.03 

10.93 

10.02 

10.0 

A15 

0.0 

.31 

-.04 

-.00 

0.0 

A16 

0.0 

-1.30 

.96 

.02 

0.0 

A21 

0.0 

-.00 

.00 

.00 

0.0 

A22 

0.0 

-2.18 

-2.00 

-2.00 

-2.0 

A23 

0.0 

-  .38 

.02 

.00 

0.0 

A24 

0.0 

.012 

-.01 

-.00 

0.0 

A25 

0.0 

1.40 

2.00 

2.00 

2.0 

A26 

0.0 

.01 

-.01 

-.00 

0.0 

A31 

0.0 

.27 

-.03 

.00 

0.0 

A32 

0.0 

-.05 

-.04 

-.00 

G.;.< 

A33 

0.0 

-1.53 

.24 

.01 

u  -  u 

A34 

0.0 

-.66 

-1.07 

-1.00 

-1.7 

A35 

0.0 

.72 

1.00 

1.00 

1.0 

A36 

0.0 

-.48 

-1.04 

-1.00 

-1.0 

A41 

0.0 

-3.23 

-11.16 

-10.02 

-10.0 
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Table  4.  (continued) 


PARAMETER 

0 

1 

4 

ITERATION  NUMBER 

10  CORRECT 

VALUE 

A42 

0.0 

-7.21 

1.38 

.02 

0.0 

A., 

43 

0.0 

29.37 

6 . 64 

9.96 

10.0 

A44 

3.0 

-7.32 

-3.67 

-4.00 

-4.0 

A45 

0.0 

.97 

-.26 

-.00 

0.0 

A46 

0.0 

-4 .10 

.49 

.00 

0.0 

AS1 

0.0 

-.05 

.05 

.00 

0.0 

A52 

0.0 

-3.46 

-5.08 

-5.00 

-5.0 

A53 

0.0 

-3 . 66 

-4.83 

-5.00 

-5.0 

AS4 

0.0 

.16 

-.03 

-.00 

0.0 

A55 

0.0 

-4.44 

-5.00 

-5.00 

-5.0 

A56 

0.0 

.13 

-.02 

-  *  CO 

0.0 

A61 

0.0 

.03 

-.04 

-.00 

0.0 

A62 

0.0 

1 . 84 

3 . 30 

3.34 

3.34 

A63 

0.0 

1.95 

3.19 

3.34 

3.34 

A64 

0.0 

-.09 

.01 

-.00 

0.0 

A6S 

0.0 

.29 

.04 

.00 

0.0 

A66 

0.0 

-6.40 

-9.75 

-10.00 

-10.0 

ail 

0.0 

3.07 

-.40 

-.00 

0.0 

812 

0.0 

-.31 

.24 

.00 

0.0 

D21 

0.0 

.00 

.00 

.00 

0.0 

022 

0.0 

.39 

.00 

-.00 

0.0 
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Table  4.  (continued) 


PARAMETER 

0 

1 

4 

ITERATION 

10 

NUMBER 

CORRECT 

VALUE 

B31 

0.0 

-.36 

.07 

.00 

0.0 

Q32 

0.0 

.24 

-.02 

.00 

0.0 

B41 

0.0 

5.85 

10.82 

10.02 

10.0 

B42 

0.0 

-.94 

.21 

.00 

0.0 

BS1 

0.0 

.06 

-.02 

-.00 

0.0 

B52 

0.0 

3.87 

5.00 

5.00 

5.0 

B61 

0.0 

-.03 

.03 

.00 

0.0 

B62 

0.0 

.06 

-.02 

-.00 

0.0 

Example  #3  -  Generalization  of  Algorithms  A  possible 


generalization  to  allow  identification  when  the  entire  state 
is  not  measured  is  indicated  by  an  example.  The  procedure 
is  to  define  a  new  state  space  in  which  the  model  parameters 
may  be  computed  by  using  output  data  and  the  output  data 
integral. 

Consider  the  second  order  phase  variable  system  in 
equations  (53)  -  (55) . 


(53)  x,  =  x, 

{54)  x2  =  a2lXi  +  a22x2  +  bu 

(55)  y  =  x1 

t  t  t 

Now  define  new  variables  w1=J  xl  (T) d-,w2=J  x2(T)dT,v=j  u(T)d 

By  integrating  (53)  -  (54)  and  substituting  w^,  w2,  the  fol¬ 
lowing  state  equations  in  the  new  states  arise. 

(56)  wi  =  w2  +  xx  (0)  -  y 

(57)  w2  =  n2lwl  n'  X22W2  +  b2V  +  x2  ^ 


The  unknown  quantities  are  a.,^,  a22,  ^2'  x2  ^  ’ 

For  the  special  case  x2 (0)  =  0,  the  model  is, 

I 

(58) 


"l\ 

J°  1  1 

iiwi| 

l 

o' 

(0) 

■  *2  > 

ia21  a22^ 

W+ 

0 

hi 

lv 

The  state  variables  of  the  above  system  are  known 
t 

since  y,  J  y(r)dT,  x^(0)  are  known.  Application  of  the 
basic  algorithm  is  now  done  straightforwardly. 
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Example  Problem  #4:  Effects  of  Model  Order  Error. 
The  robustness  of  the  model  adaptive  algorithm  to  model 
order  errors  is  indicated  by  the  present  example.  A  sec¬ 
ond  order  model  is  fitted  to  data  generated  by  a  third 
order  system  using  a  slight  modification  of  the  basic 
algorithm.  The  modification  is  that  the  identification 
is  not  stopped  when  an  A,B  matrix  is  determined.  That  is, 
the  basic  algorithm  is  applied  to  e.  data  matrix  (XO,  UO) 
using  n+p  measured  state  and  input  vectors  and  an  a»B  is 
determined  which  fits  th$  data.  Every  vector  is  then  dis¬ 
placed  to  the  left,  and  a  new  data  vector  replaces  the 
last  column  of  the  matrices  XO,  UO.  The  algorithm  is 
applied  to  the  new  data  matrix.  Previously  determined 
model  matrices  A,B  are  used  as  starting  points. 

The  third  order  system  chosen  is  a  model  of  a  gun 
turret  for  use  on  ships.  It  is  described  by  the  transfer 
function  shown  below. 

(  c  q  v  _  _ 1 2  ♦  1 

'  '  V  "  S  (S+18XS+3.5) 

0  is  the  turret  position  and  V  is  the  voltage  applied  to 
the  turret  drive.  A  step  of  three  volts  was  used  as  the 
drive  input. 
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Table  5  shows  that  the 
in  the  sense  that  it  converges 
lower  order  model  is  fitted  to 
order  system.  For  the  example 
are  close  to  those  obtained  by 
damped  pole. 

Table  5.  Second  Order  Model 
Iterations 

Fit  to  a  a 

No.  Convergence  all  a12 

1  3  -.001  1.006 

2  3  -.002  .994 

Coefficients  obtained 
by  disregarding  pole 
at  -18 

0.0  1.000 


basic  algorithm  is  robust 
to  a  unique  answer  when  a 
data  generated  by  a  higher 
chosen,  the  coefficients 
neglecting  the  most  heavily 

Fitted  to  Third  Order  Data 

a21  a22  bl  b2 

-1.182  -2.61  -.000  .532 

-.151  -3.243  .001  .575 

0.0  -3.5  0.0  .535 
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SUMMARY 


A  model  adaptive  identification  scheme  was  developed 
which  yields  the  system  and  input  matrices  of  linear,  time 
invariant  multi-input  systems.  Measurements  of  the  state 
and  inputs  to  the  system  were  assumed. 

Solution  of  linear  equations,  generation  of  the 
model  transition  matrix  and  its  integral,  and  integration 
of  the  model  equations  were  the  primary  numerical  tech¬ 
niques  required. 

The  primary  requirement  for  unique  convergence  was 
the  non-singularity  of  the  data  matrices. 

Some  rules  of  thumb  on  choice  of  sampling  interval 
were  given. 

Several  example  problems  showed  that  the  algorithm 
is  rapidly  convergent. 
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CHAPTER  lV 


PARAMETER  ESTIMATION  BY  NON-LINEAR  OBSERVERS 

INTRODUCTION 

The  previous  chapter  assumed  that  all  states  were 
measured  and  developed  a  model  adaptive  identification 
scheme  which  worked  very  well  for  that  case.  In  many  cases 
of  engineering  interest,  it  may  indeed  be  possible  to  meas¬ 
ure  all  these  states.  However,  if  the  states  are  not  all 
measured,  then  a  technique  is  desired  for  parameter  estima¬ 
tion  which  will  work  using  the  measured  states  which  are  a 
subset  of  the  system  states.  The  present  chapter  discusses 
parameter  estimation  using  a  generalization  of  linear  ob¬ 
servers.  Preliminary  background  material  is  developed  for 
linear  observers  and  the  effects  of  unknown  parameters  are 
discussed.  Following  the  background  discussion,  the  theory 
of  non-linear  observers/parameter  estimators  is  developed. 
Finally,  a  numerical  example  is  presented. 

LINEAR  OBSERVERS 

At  the  present  time,  linear  observer  techniques  used 
for  estimating  the  states  of  linear  systems  are  Well  known, 
principally  through  the  work  of  Luenberger  (SE-2] .  In  gen¬ 
eral,  observers  estimate  states  with  an  error  which  asymp¬ 
totically  approaches  zero.  The  design  of  these  observers 
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is  predicated  upon  the  assumption  that  one  knows  the  system 
and  the  input  matrices  for  the  system  to  be  observed.  Some 
things  about  linear  observers  which  are  not  well  known,  how¬ 
ever,  are  the  effects  of  system  and  input  matrix  changes  on 
the  state  estimate.  The  following  brief  development  shows 
the  effects  of  parametric  changes  in  the  observed  system  on 
the  observer  states.  The  result  of  the  development  shows 
that  the  error  in  a  state  estimate  is  dependent  not  only  on 
the  observed  system's  state  and  input  to  the  system,  but 
also  on  the  variations  in  the  system  and  input  matrix.  It 
is  assumed  that  the  observed  system  is  described  by  equa¬ 
tions  (60)  and  (61)  where  C  is  the  known  output  matrix. 

(60)  x.  “  A  x  +  Bvv,  x(G)  not  known 

(61)  X  -  C  x 


Further,  assume  that  the  system  is  erroneously  thought 

to  be  modeled  by  equations  (62)  and  (63) .  Consequently,  the 

observer  will  be  designed  using  the  erroneous  system  matrix 

h  and  the  erroneous  input  matrix  Q  . 
m  m 


(62) 

(63) 


x  =  A  x  + 
— tn  m— m 


D  vi, 
m— 


x  (0)  not  known 
— m 


X 


Cx 


Equation  (64)  describes  the  observer. 


(64)  w  *  Fw  +  G£  +  Hu 
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The  desired  objective  is  to  have  the  observer  state 
w  =  Tx  +  e  where  _e  goes  to  zero  at  some  satisfactory  rate. 
For  the  case  at  hand,  T  is  a  non-singular  transformation 
and  we  would  like  to  take  x  =  T~^  (w)  as  the  state  estimate. 
If  che  observer  is  designed  properly,  the  error  will  in  fact 
approach  zero  asymptotically.  By  performing  some  algebraic 
manipulations  and  subtracting  equation  (60)  from  equation 

(64)  the  expression  for  the  error  in  the  state  estimate  is 
given  by  equation  (65). 

(65)  T-1e  »  T_1FT  T_1e  +  (T-1FT  +  T^GC-A^AA)  x 

+  (T“1H-B  -  AB) u 
m  ~ 

If  the  state  estimate  is  to  asymptotically  decay  to 
zero,  we  need  to  eliminate  the  forcing  terms  and  further, 
ensure  that  the  observer  is  stable.  Since  it  is  thought 
that  equation  (62)  models  the  system,  we  oncost  the  matrices 
T,  F,  G,  H  to  satisfy  equations  (66*  and  (67). 

(66*  TA  -  FT  »  GC 
tn 

(67)  H  -  TBm  a  0 

Choosing  the  observer  to  satisfy  those  equations, 
yields  the  error  in  the  state  estimate  shown  by  equation  (68) 

(68)  T-1e  =<  (AW-T_1GC)  T-1e  -  A  Ax  -ABu 
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It  is  obvious  that  the  state  estimate  error  is  com¬ 
posed  of  two  parts.  The  first  part  is  the  free  response 
resulting  from  the  unknown  initial  condition  x(0)  since,  if 
it  were  known,  we  could  set  £(0)  =  x(0)  -T"1w(0)  =  JD  and 
have  no  free  response  term.  The  second  is  the  forced  response 
caused  by  the  unknown  parameter  variations  AA  and  Ab.  These 
parameter  variations  enter  in  a  multiplicative  fashion,  and 
their  effect  on  the  state  estimate  depends  entirely  on  how 
the  system  is  excited  and  what  these  parameter  variations 
are.  Utilization  of  the  observer  to  provide  the  state  esti¬ 
mate  for  the  model  adaptive  identification  scheme  of  the 
previous  chapter  will  provide  a  state  estimate  which  is  in 
error  by  the  amount  shown  in  equation  (68)*  Consequently, 
the  linear  observer  has  not  been  found  to  be  useful  for 
identification.  Modification  of  the  way  in  which  the  system 
equation  is  written  and  considered  yields  a  non-linear  ob¬ 
server  whose  state  vector  is  the  original  system  state  and 
the  parameter  vector. 

NON-LINEAR  OBSERVERS 

The  last  section  developed  the  expression  for  the 
error  in  the  state  estimate  of  an  observer  improperly  de¬ 
signed.  If  the  observer  for  a  linear  system  is  designed 
based  on  erroneous  system  and  input  matrix  information,  the 
observer  estimates  the  state  of  the  observed  syst  ..  witv  an 
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error  which  is  defined  by  linear  differential  equation.  This 
error  is  uncorrectable  since  it  depends  on  knowledge  of  the 
unknown  matrices  AA,  AB,  and  the  state  which  is  also  unknown. 
A  method  of  circumventing  the  difficulty  of  unknown  system 
parameters  and  system  states  has  been  developed  using  non¬ 
linear  observers.  Essentially,  what  is  done  is  to  rewrite 
the  linear  system  of  equations  in  non-linear  form  with  an 
augmented  state  vector.  This  new  non-linear  system  has  as 
its  states  the  original  state  x  plus  a  vector  of  parameters 
for  the  system  matrix  A  and  for  the  input  matrix  B.  The 
parameter  vector  for  the  system  matrix  is  called  and  it 
is  formed  by  stacking  the  rows  of  the  system  matrix  A  into 
a  vertical  vector.  In  similar  fashion,  the  parameter  vec¬ 
tor  b  for  the  input  matrix  B  is  formed  by  stacking  the  rows 
of  the  B  matrix  into  a  vector.  Using  the  Kronecker  product 
the  non-linear  system  then  is  described  by  equations  (69) 
through  (72) .  Notice  that  constraints  on  the  parameters 


are 

also  included.  The  constraints  are 

the  e .  . ,  d 

ij 

h.  . , 

13 

terms. 

(69) 

x  =  I  Q  +  I  Q  uS) 

(70) 

a  =  0  e  ,  <a .  .  <d .  .  i 

-  -  13  13 

,j  -  1,2,. 

(71) 

b  =  0  k,  . <b .  .  <h .  . 

-  -  13  13  13 

i  0  1,2,. 

(72) 

Z  =  Cx 

j  -  1,2,. 
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Given  the  non-linear  augmented  system  we  construct  an 
observer  which  will  yield  the  entire  augmented  state,  that 
is,  the  original  system's  state  plus  the  parameter  vectors 
a,  b.  To  construct  the  observer,  we  build  in  a  structure 
similar  to  that  in  the  augmented  system  and  include  parame¬ 
ter  constraints  if  they  are  known.  The  general  form  of  the 
observer  is  shown  in  equations  (73)  through  (75) . 

(73)  Hj_,cF^(t)^  +  +  Ia— ^ 3  +  Gi 

(74>  w2=F22(t)-2  +  F2i{t*-1  +  G2{t^ 

(75)  ”3  =  F33(t)w3  +  F3l(t)wx  +  G3(t)£ 

It  is  desired  that  the  observer  state  vector  corres¬ 
pond  to  the  augmented  state  vector  of  the  original  system 
within  some  error.  The  notation  involved  is  shown  in  equa¬ 
tion  (76)  through  (78)  • 

(76)  w3  a  x  +  e_3 

(77)  w2  a  £  +  e_2 

(78)  w3  =  b  +  e3 

As  in  the  last  section,  it  is  desired  to  develop  the 
error  equations  for  the  non-linear  observer.  After  some 
algebraic  manipulation  and  subtraction  of  equations  (69) 
through  >72)  from  equations  (73)*  (74)  and  (75),  it  is  seen 
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that  the  error  equation  is  given  by  equations  (’jpg)  through 
(81)  . 

(79)  e ^  =  (Fj^  +  i  +  +  +  3  + 

(Pll+  GjCjx 

(QO)  e2  =  ^22— 2  +  F21-l  +  F22—  +  ^F2i  +  G2C^- 

(81)  £3  =  F33£3  +  ^3^2.  +  F33^  +  (F3i  +  F3C)* 

Examining  equations  (79)  through  (81)  it  is  seen  that 
the  error  equations  are  driven  by  the  state  which  is  unknown. 
To  eliminate  the  state  from  the  error  equations,  we  intro¬ 
duce  the  constraints  shown  by  (82)  through  (84). 


(82) 

PU  = 

-GiC 

(83) 

F21 

G2C 

(84) 

F31  = 

-GjC 

The  £in?.l  form  of  the  error  equation  is  shown  in 
equations  (85)  through  (87). 


(85) 

k  = 

(A  -G1(t)C)e1  +  IQw^ lc2  +  ISn^e 

(66) 

k  D 

E22(t)e2  -Gjtt)^  ♦  F22(t)a 

(87) 

i3  - 

F33(t)«3  -G3(t)Ce1  +  F33(t)b 
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In  general,  it  may  be  highly  desirable  to  let  the  gain 


matrices  G^,  G2,  be  time  varying.  Therefore,  the  feed¬ 
back  matrices  F21'  F31  w-*-H  also  b®  time  varying.  If 

the  observer  is  to  act  as  a  successful  state  estimator  and 
parameter  estimator,  the  error  equation  must  eventually  go 
to  zero.  In  general,  this  means  that  the  error  equation 
must  be  stable  and  the  matrices  F22(t),  F^ft)  must  go  to 
zero  in  such  a  fashion  that  the  stability  of  the  error  equa¬ 
tions  is  preserved.  The  matrices  F22,  must  go  to  zero 
because  inspection  of  (86)  and  (87)  indicates  that  e^,  e^ 


are  forced  by  some  linear  combination  of  the  parameter  vec¬ 
tors  a_,  b.  These  forcing  terms  must  eventually  be  elimi¬ 
nated,  and  since  _a,  b  are  constant,  the  forcing  terms  can 
only  be  eliminated  by  letting  F22,  F ^  eventually  approach 
zero  as  the  state  and  parameter  estimates  converge  to  their 
correct  values.  In  general,  the  observer  is  non-linear  and 
time  varying,  hence,  its  stability  properties  can  be  criti¬ 
cally  dependent  on  the  input  u  and  the  state  of  the  observed 
system  (w^  -  +  £^)  thus,  the  scheme  may  work  for  some 

selection  of  inputs  and  not  for  others.  As  is  common  for 
other  identification  schemes,  all  state  variables  must  be 
excited.  Another  consideration  is  that  the  stability  of 
the  observer  cannot  be  assessed  prior  to  the  computation. 

The  cause  of  the  problem  is  that  the  error  equations  contain 
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the  unknown  matrices  A,  B,  as  well  as  the  unknown  state. 


Selection  of  the  gain  matrices  G- ,  G2*  G^  must  be 
done  with  some  care.  In  contrast  to  the  linear  observer, 
where  elements  of  the  gain  matrix  would  be  non-negative, 
the  gain  matrix  elements  are  either  positive  or  negative. 

A  study  of  the  model  structure  of  the  identified  system  is 
required  to  determine  the  proper  sign  on  the  gain  elements 
for  the  parameter  vector  estimate.  For  example,  suppose  a 
first  order  system  with  two  parameters  a,  b  is  to  be  identi- 


fied. 

The  system  and 

observer  are  described  in 

(88) 

x  =  -ax  +  bu 

a,  b  >  0 

(89) 

a  =  b  =  0 

(90) 

y  =  x 

(91) 

Wf  =  w2wx  +  w3u 

+  g^  (y-w1)  w2,  w3  >0 

(32  . 

w2  =  9 2 (y 

(93) 

w3  =  93(y-w1) 

The  gain  g^  should  be  positive  to  insure  that  has 
negative  feedback.  Gain  should  be  negative.  The  reason 
for  this  is  that  if  y-w^  is  positive,  then,  x  in  greater 
than  .  Now  assume,  temporarily,  that  £  b.  If  x  > 

then  w2  is  too  large  and  should  be  decreased,  i . c . , <  0, 
or  q,  <0.  If  v -w,  <  0,  then,  x  <  w,  .  when  x  <  w  ,  then 
for  w  -  b,  w.,  is  too  small  and  should  lv-  inoreas  -n ,  i.e., 

J  *»  ’ 

>  0  or,  again,  g2  <0. 
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By  similar  reasoning,  it  is  clear  that  g3  should  be 
positive. 

The  determination  of  the  gain  sign  for  higher  order 
systems  is  a  generalization  of  the  above  example.  Each 
term  is  considered  by  itself,  i.e. ,  all  parameters  but  one 
are  assumed  correct  in  order  to  determine  the  relationship 
of  the  sign  of  y-w^  and  the  parameter  gain  sign  to  achieve 
the  required  sign  on  the  parameter  estimator  rate. 

OBSERVER  STABILIZATION 

Stabilization  of  the  observer  is  difficult  because 
constant  feedback  tends  to  make  the  parameter  estimates 
diverge  from  their  true  values.  Such  divergence  causes 
observer  output  error  which  in  turn  causes  the  parameter 
estimates  to  be  in  error  and  the  cycle  tends  to  repeat.  To 
illustrate  the  point,  the  error  equations  associated  with 
the  previous  example  are  given  below.  A  parameter  feedback 
term  is  added  in  (92)*  (93). 

(94)  -•  •  (g  ^  +  v'2^c'i  *  °tx  +  C3U 

(9?)  e,  *  f->c'2  +  ^?2C1  +  ^2a 

(96)  e.  a  -f,c  +  go,  +  £.,b 
J  3  3  3  1  j 
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Analysis  of  (95)  and  (96)  indicates  that  e0  and  e, 

tend  to  approach  a  and  b  if  f ^  and  f^  are  positive  constants 

If  £2  ~  £3  =  then,  a  equilibrium  solution  to  (94)  -  (96) 

is  e,  =  e0  =  e-.  =  0.  The  difficulty  is  the  solution  is  not 

stable.  A  choice  of  f^  and  f_,  must  be  made  such  that  f_ 

Z  o  2 

and  f-j  tend  to  zero  while  maintaining  the  system  stability. 

A  stabilization  technique  which  has  worked  is  to  modi¬ 
fy  the  parameter  equation  so  that  the  parameter  rate  is  zero 
at  the  end  of  some  time  interval  during  which  the  observer 
has  provided  acceptably  accurate  parameter  estimates.  The  • 
parameters  are  denoced  by  p^,  i  =  l,2,...,n  where  n  is  the 
total  number  of  parameters.  One  possible  parameter  equa¬ 
tion  is  (97 )  . 

(97)  pt  =  -  a  e"t'Tpi,  pi(o)  =  Pio  i  =  l,2,...,n 
The  time  solution  for  p^  is  (98) 

(98)  pi  =  Pio  EXp(a  t  (EXP(-t/T  -1));  i  =  l,2,...,n 

There  are  two  adjustable  parameters  which  relate  the 

steady-state  value,  p^(<°)  to  the  initial  value  P^Q.  The  t 

« 

parameter  primarily  controls  the  rate  at  which  p^  approaches 
zero,  and  the  a  parameter  primarily  controls  the  ratio 
Pj_(ro)  /Pj_0*  The  exact  expression  for  p^(“)/p^Q  is  (99). 
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p,  (•) 

(99)  -  =  EXP(-<xt)  i  =  1,2,  ...  ,  n 

pio 

Let  P^o  be  the  correct  value  of  the  parameter,  then, 

p^(“)  is  the  approximate  value  of  the  parameter. 

P±  <00) 

If  <Xt  =  .02,  then,  - -  =  .98,  and  the  parameters 

l 

of  the  approximate  model  are  only  two  percentage  points  dif¬ 
ferent  from  the  parameters  of  the  exact  model.  Application 
of  the  approximate  parameter  approach  to  the  previous 
example  leads  to  (100)  -  (101)  as  estimators  for  the  a,b 
parameters. 

(100)  w2  «  -ae-t//Tw2  +  g^y-w^ 

(101)  w3  =  -ae-t//Tw3  +  g3(y-w1) 

The  estimators  are  stable  for  a  >  0  during  the  time 
period  for  which  the  feedback  is  effective.  By  proper  de¬ 
sign  choice  of  t,  g^,  g2,  g3,  the  observer  can  be  made  to 
give  a  stable,  accurate  estimate  of  the  system  parameters. 


NUMERICAL  EXAMPLE 


A  numerical  example  to  illustrate  the  technique  was 
worked.  Equations  (102)  “  (ll4)  describe  the  second  order 
system  and  its  fifth  order  state  and  parameter  estimator. 

(102)  -  x2 

(103)  x2  =  a21xL  -  a22x2  +  bu 

(104)  y  = 

(105)  a21  -  a22  =  b  =  0 

(106)  a21  =  4 

007)  a22  «  2 

(108)  b  «  1 

(10S)  w1  »•  w 2  +  (y-w^) 

(110)  w2  «  "w3wl  ‘W4W2  +  W5U  +  g2 

(111)  «  -pw.j  -  400  (y-Wj) 

(112)  w.  -  -}iw4  -  1000  (y-wx) 

(113)  a  -pV*.  +  ^0  ^y-w|) 

(114)  P  a  .4  EXP(-20t) 

The  state-parameter  estimator  to  system  correspondence 
was  v/^ ,  %■',2'*x2,  wns'*a2i»  w4"*a2'?'  Table  6  shows  the 

numerical  results . 
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Table  6  .  Observer  Identification  Example 


Time 

(sac) 

Y"W! 

W3 

W4 

W5 

0 

-4. 22xl0~3 

2.50 

3.50 

1.50 

.4 

-5.43xl0~4 

2.36 

4.10 

1.19 

.8 

-7. 14xl0~4 

2.15 

3.96 

1.09 

1.2 

-7.31xl0~4 

2.14 

4.01 

1.06 

1.6 

-4.86xl0-4 

2.13 

4.06 

1.03 

2.0 

-2 . 02xl0~5 

2.11 

4.05 

1.02 

2.4 

-1 . 85xl0-5 

2.12 

4.06 

1.02 

2.8 

-6 . 16xl0-5 

2.11 

4.07 

1.01 

3.2 

7 . 28xl0~5 

2.09 

4.01 

1.01 

3.6 

-1 . 3.1xl0~4 

2.08 

4.01 

1.01 

4.0 

-2 . 83x10 

2.08 

4.02 

1.00 

4.4 

3 . 35x10  5 

2.06 

4.00 

1.00 

4.8 

-1. 20xl0-4 

2 . 06 

3.99 

1.00 

5.2  ' 

-1. 58x10 “5 

2.06 

3.99 

1.00 

5.6 

1.76x10'“  5 

2.05 

3.95 

1.00 

5.8 

-9. 93xlQ~5 

2.04 

3.96 

.99 

Correct  Value 

0 

2.00 

4.00 

1.00 
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The  estimator  converged  to  the  correct  values  rapidly 
for  the  choice  of  initial  conditions  and  parameter  stabili¬ 
zation  technique.  It  is  felt  that  the  estimation  scheme 
is  well  suited  to  an  analog  mechanization. 


SUMMARY 

A  combined  state-parameter  estimation  technique  was 
developed.  The  scheme  is  based  on  a  generalization  of  the 
linear  observer  concept  and  was  shown  to  converge  well  for 
a  specific  example. 

The  estimator  is  multi-linear  and  can  be  unstable  un¬ 
less  a  time  varying  feedback  on  the  parameter  estimation 
portion  is  used.  A  method  of  stabilizing  the  observer  was 
presented. 

Computationally,  the  estimator  is  good  because  it  can 
be  implemented  sVraightforwardly  on  an  analog  computer. 
Feedback  gain  adjustment  is  relatively  easy.  Since  no  ana- 
log-to-digital  converters  are  required,  on-line  simulation 
should  be  a  potential  use  of  the  technique. 
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CHAPTER  V 


ESTIMATION  OP  PARAMETER  SUBMATRICES 

INTRODUCTION 

The  present  chapter  is  concerned  with  determining  the 
parameter  submatrices  of  network-like  physical  systems. 

These  parameter  submatrices  are  contained  in  the  system  and 
input  matrices  of  the  system  which  are  assumed  to  be  known. 
A  further  consideration  is  determining  if  enough  measure¬ 
ments  have  been  taken  to  uniquely  define  the  parameter  ele¬ 
ment  matrices. 

As  is  well  known  [Koenig,  FE-8]  if  there  are  e-ele¬ 
ments  in  the  system  and  if  two  e-measurements,  i . e. , 
e-current  measurements  and  e  -voltage  measurements  in  the 
case  of  electrical  networks,  are  made,  then  each  pair  of 
measurements  defines  a  parameter  such  as  a  resistor  or  a 
capacitor.  By  using  the  generalization  of  Kirchoff's  volt¬ 
age  laws  systematically  through  the  theory  of  linear  graphs, 
a  sufficient  number  of  measurements  to  estimate  the  values 
of  the  e-parameters  is  e.  The  techniques  of  the  present 
chapter  allow  further  reductions  in  the  number  of  measure¬ 
ments  by  making  full  use  of  the  known  topological  informa¬ 
tion  concerning  the  system  and  the  results  from  the  identi¬ 
fication  methods  previously  presented. 
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The  state  equations  are  assumed  to  be  in  the  form  as 
given  in  the  Appendix  and  it  is  further  assumed  that  the 
network  does  not  have  excess  dynamic  elements,  e.g.,  capaci¬ 
tors  in  parallel  and  inductors  in  series.  Mutual  inductance 
is  allowed  under  the  last  assumption.  It  is  also  assumed 
that  the  A,B  matrices  have  been  determined  by  some  identi¬ 
fication  technique  and  that  there  are  no  non-state  connected 

sources,  i.e.,  K  K  .  are  null  matrices. 

ra  ga 

Under  the  above  assumptions,  the  known  information  is 
summarized  in  the  following  equations. 

(115)  x  =  Ax  +  Bu 

Where  x(t) ,  u(t) ,  tc[0,T]  are  known,  and  A,B  are 


known. 

(116) 

,J11  0  |/A11A12\ 

/FU  F12 

\°  J22/  \A21  A22 

iF21  F22 

/  Jn  0  \  /B1\  = 

'Hl\ 

(117) 

\°  J22 1  '  B2 ' 

>h2  ' 

(118) 

(Jn)  =  (cs) 

(119) 

(J22)  - 

(120) 

Pu  =  -T2(R+tJ  t/ 

(121) 

f12  -  t2<r+t4  a-1**)-1  t*  g- 

'1T3  +  T. 

(122) 

F21  =  A 

-  T^G+^R”1^)"1 

T4R"1t2 

(123) 

F22  =  -T3 

(g+t4r_1tJ) _1  t3 

(124) 

H1  =  Kcd 

(125) 

H2  *  Kcd 

The  unknowns  in  (116)  -  (125)  are  C  ,  L  ,  R,  G,  but  it 

s  s 

is  known  that  Cg,  R,  G  are  diagonal,  and  Lg  is  symmetric  and 
positive  definite. 

SOME  BASIC  RESULTS 

The  first  results  concern  the  independence  of  the 
matrix  equations  (120)  -  (123) . 

Result  #1  F12  -  F21t  =  2TX 

To  see  this,  use  the  matrix  identity  ABfc (C+BAB^) = 
(A-1  +  B^^B)"1  bV*1  and  rewrite  (122)  as, 

(126)  F21  =  -T1  -  T2 (RtT^G_1T4) '1tJ  G-1^ 

•  fc 

and  subtracting  F21  from  F^7,  the  result  follov;s. 
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Result  #2 


F12+F21  =  2T2(R+tJ  G'1^)  h^G'1 

The  result  follows  from  above. 

The  next  results  are  concerned  with  the  relationship 
of  the  number  of  resistive  and  dynamic  elements  to  the 
singularity  of  the  topology  matrices  and  may  be  established 
by  simple  dimension  analysis  of  these  matrices. 

Result  #3  det  T^  ^  0  only  if  nc  ®  n^ 

Result  #4  det  T-  ^  0  only  if  n  =  n 

Result  #5  det  T,  ^  0  only  if  nc  =  n 

<3  A  ^ 

Result  #6  det  T.  /  0  only  if  n  »  n 

*#  r  g 

Since  the  matrices  R,G,  are  diagonal,  the  results 
above  imply  the  following  results. 

Result  #7  1'2  is  non-singular  only  if  the  number  of 

tree  capacitors  is  equal  to  the  number  of  chord  resistors. 

Result  #8  is  non -singular  only  if  the  number  of 

tree  conductors  is  equal  to  the  number  of  chord  inductors. 
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DETERMINATION  OF  C  ,  L 

3  _ S_ 

In  general,  the  dynamic  element  matrices  will  not  be 
determined  by  state  and  input  measurements  only.  Some  ad¬ 
ditional  measurements  may  be  required,  hence,  the  purpose 
of  the  present  section  is  to  develop  a  technique  for  using 
the  information  available  from  the  state  and  input  measure 
ments  to  full  advantage  and  to  show  how  to  select  the  addi 

tional  measurements  required  to  uniquely  define  C  ,  L  . 

s  s 

The  off  diagonal  blocks  Aj^,  A22  of  the  system 
matrix  A,  and  the  input  matrices  B^ ,  B2  are  used  to  de¬ 
termine  as  many  elements  of  C  ,  L  as  possible,  and  the 

s  s 

remaining  elements  must  be  specified  by  a  combination  of 
topological  information  and  additional  measurements. 

Using  the  constraint  between  F^*  F2^  and  noting 

that  F^2  =  csA12  anci  ^21  a  ^-^21'  fc^e  tlirPe  m^trix  equa¬ 
tions  available  for  determining  C  ,  L  are, 

5  8 


(127) 

(128) 
(129) 


(131) 


CsA12  ~A21Ls  "  2T1 


C  B,  -  H. 
S  1  1 


LsB2  =  H2 

Two  sets  of  equations  may  be  derived  from  (130)  - 


by  considering  either  L  or  C,  as  independent. 

s  s 
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t.-':  tv.  . 


and 


U30)  <VA12B2:V  =  (2T1B2  +  A21H2:H1> 


U31)  VA2iBiiB2>  -  <2tJbx  + 


Now  drawing  on  the  diagonal  structure  of  Cg 
letting  X  =  (A^BjjBj.),  Y  =  +  A  ,  (130)  can 

be  written  as  (132)  . 


(132)  C.x.  y.  ■  ,  0 

i— 1*  =  i.1*  i  =  1,2,  ...  ,  n 


By  adding  the  components  of  x.  ,  y,  ,  the  C.'s  can 
be  found  from  (133). 


n  +  P 
l  *c 


*»1  +  P, 


(133)  Ci  £  j  =  £  Yij  i  =  1,2,  ...  ,  nc 

5-1  5-1 

Let  r  be  the  number  of  non-zero  row  sums  of  X,  then, 

i  ~  1,2,  ...  ,  r 

(134) 


nf  - 

pc 

r 

1  =  1 

nl  * 

Pc 

r 

j  =  l 

i5 

The  number  of  elements  of  C  which  must  be  specified 

s 

by  additional  measurements  is  n-r. 


b8 


Since  L  may  be  non-diagonal,  the  determination  of  the 
s 

elements  of  its  elements  is  more  difficult. 

The  Kronecker  product  can  be  used  to  advantage,  how¬ 
ever,  by  post-multiplying  (130)  by  (A21B1jB2)t  to  convert  to 
a  system  of  square  matrix  equations. 

Let  W  =  (A21Bi:b2)  (A21Bi:B2)t  and  Z  «  (-2T1tB1  + 

H2)  (A21B1*B2^ then'  ^129^  is  written  in  the  Kronecker  pro¬ 
duct  notation  as  (135) . 


(135)  I«W  t  =  z  * 


0 


*w 


is 


2  i 

where  W  is  a  n.xn  and  i  is  an  n,  -  dimensional  vector  cora- 
l  —s  l 

posed  of  the  transposed  rows  of  Ls» 

If  has  rank  n^,  then  the  problem  is  solved,  and  the 
elements  of  Ls  are  given  by, 

(136)  e  =  W-t®i  z 

The  more  usual  case  is  that  rank  w  »  r  <  n  .  and  since 

x 

h 

rank  Af®  =  rank  A'*  rank  B,  rank  (I*v)  »  rank  (I«w  )  =•  n  .r 
2 

and  -  r*n  components  of  must  be  specified  by  topo¬ 
logical  information  and  additional  measurements. 

Now  drawing  on  the  symmetric  structure  of  L  ,  there 

S 

arc  n/.  i  +  ^  possible  different  elements  of  t_.  Define 
2  ** 
the  reduced  matrix  and  vector  D,  z  of  order  n  ^  x 

(n  +1)  and  n  t  ^ n 4^ \ |  x  1  respectively  obtained  from 
l  2 
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by  deleting  redundant  entries  caused  by  the  symmetry  of 
Ls  by  (137)  . 

(137)  (D)^  =  z 

The  equation  can  be  further  reduced  in  order  by  per¬ 
muting  rows,  columns,  and  components  of  D,  to  jz  to 

yield  the  following  equation. 

(138)  (Dr*r*;Dr*25  i-fi-l  =  -1 

-2 

n,(n  ,+i) 

where  r*  =  rank  D  and  D  is  an  r*x  - 2 -  matrix.  Solving 

(138)  for  j,y 

(139)  ^  a  iL]_  ~  Dr*r*  -&z 

hence,  the  vector  must  be  specified  either  by  topologi¬ 
cal  information  or  additional  measurements  in  order  to  de¬ 
fine  the  complete  vector  uniquely.  The  number  of  these 

>  *.  n,(n,+l) 

measurements  rs  i  j _ -r. 

2 

DETERMINATION  OF  Rf  G  SPECIAL  CASE  T 4  -  (0) 

For  the  special  case  where  there  is  no  tree  to  cotree 
resistive  coupling,  the  determination  of  R,  G  is  greatly 
simplified. 

It  is  assumed  that  C  ,  L„  have  been  determined, hence, 

s  s 

from  that  knowledge  and  knowledge  of  the  system  matrix  A, 
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the  submatrices  F^  -  CsAll'  F22  ~  LSA22  ma^  determined> 
Then  using  the  topological  knowledge  of  the  system,  R,  G 
are  defined  in  (140)  and  (141) . 


(140)  Fn  *  -T2R~1T2t 

(141)  F22  =  -TitG~1T3 

2  2 

Converting  (140)  and  (141)  to  n  and  n  dimensional 

^  y 

vector  equations  by  the  Kronecker  product  after  pre  and 

£ 

post-multiplying  by  T2>  as  appropriate, 

(142)  -T2tT2^i’9tT2x  »  £* 

(14J)  -T^TjSTjTj^  -  fj2 

where  x,  y  are  composed  of  vectors  of  the  rows  of  R~^,  C--i, 
respectively,  and  .£*22  are  vectols  of  the  rows  of 

T2tFiiT2.  TjF22Tjt'  respectively. 


The  number  of  components  of  x,  y  which  must  be  speci- 

2  2  2  2 

fied  to  define  x,  y  uniquely  is  n  -  r~  ,  n  -  r  where 
—  ■*-  i  i  g  g 

r2‘  r3  are  t*ie  ranks  of  T2,  respectively. 


APPROXIMATE  SOLUTIONS 

It  may  happen  that  one  is  interested  in  the  approximate 
values  of  R,  G,  Cg,  Lg  to  serve  as  either  a  .starting  point 
for  more  exact  solution*,  as  in  the  case  of  successive 
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approximations,  or  because  approximate  values  are  sufficient. 
The  generalized  in/erse  may  be  used  to  great  advantage  be¬ 
cause  of  the  special  structure  of  the  topology  matrices, 
i.e.,  their  elements  are  0,  ±  1  and  the  construction  of  che 
generalized  inverse  is  much  easier  than  for  the  general 
case. 


By  using  the  matrix  identity  (C+S^A  "^B)  ^  ~ 
(A+BC~^B^) (120,  (123)  may  be  rewritten  as, 

(144)  Fn  =  -I*2 R"1^  +  T2R_1tJ  (G+T4R‘'1T4t)  "1T4R"1V2t 


t„-l*  .  *  t„-lm  /t%im  t,,-l_  ,-lm  t„-l„ 


(*45)  F22  a  -T3  G  T3  +  T3  G  T4(R+T4  G  T4)  T4  G  T- 


Let  £12  -  P12  +  F21 


3  3 

t 


(146) 


12 


=  -2T2R~1T4t(G  +  T4R-1T4) “1T3 


Therefore,  using  the  generalized  inverse, 

(147)  (G~V  =  -(I  +  T2IF12T3IT4)IT2IF11(T2t)1 

T 

(148)  (R”1)1  =  (I+T'2F12T3IT2F11(T2t)1 

*T 


Using  the  diagonal  structure  of  R, 

diagonal  elements  of  (146),  (147)  are  — 

^i 


G,  the  i-th 

,  respectively. 
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MAIN  SOLUTION  TECHNIQUE  FOR  R, G 


Since  the  identification  technique  provides  more  in¬ 
formation  than  the  A,B  matrices,  this  information  is  used 
to  help  determine  R,G.  The  outline  of  the  scheme  is  to  de 

termine  C  ,  L  by  the  techniques  presented,  and  then  use 
s  s 

the  state  equations  plus  knowledge  of  the  system  topology 
to  determine  R,G .  The  remaining  information  used  is  the 
state  derivative  which  is  easily  computed  from  x  =  Ax  +  Bu 
where  A,  x,  P,  u  are  known.  A  prime  benefit  of  this  tech¬ 
nique  is  that  linear  algebraic  solutions  are  required  in¬ 
stead  of  non-linear  solutions. 

The  technique  is  formulated  for  the  case  where  the 
independent  variables  are  the  state  variables  and  the  tree 
conductor  voltages.  The  corresponding  equations  are. 


(149)  (1’  R-1)  (-T~v  -  T. fcv  )  =  T-i  -  C  v  +  H  ,i  . 

'  2  2— s  4  -g'  x— s  s— s  cd— cd 


(150)  -T-^v  =  L  i  •  -  T.  fcv  +  K  „  ,e,  , 

'  3  -tg  s— s  1  — o  £d-od 


The  unknowns  in 

rank  T-, ,  then  n„  -  r.. 
3  g  3 

of  v  are  required  to 
~9 

then  that  v  has  been 
“9 

and  computation  d£  x, 
(151)  Vs=(vs<tc), 


(14y) ,  (150)  are  v  ,  R~^.  If  r.  = 
additional  measurements  of  components 
define  v  (t)  te(0,T]  uniquely.  Assume 
defined  and  from  the  chart  records 
form  the  following  matrices. 

...  .  ^(t^)) 


(152) 

I  = 
s 

•  is(tk>) 

(153) 

^  = 

(.  ,  , 
^(t0), 

(154) 

U  = 
s 

(icd(to^'  ** 

•  '  ±cd(V  ) 

(155) 

V  = 
9 

(VV'  •** 

•  V1*1) 

Using  (153)  -  (155) ,  equation  (148)  can  be  expanded 
to  yield  the  matrix  equation  (156) . 

(156)  T0R'“1(-T0tVo  -  T.Sf  )  =  T,  I  -  C  Y  +  H  ,U 

2  '2s  4  g'  Is  ss  cds 

The  matrix  T-Sf  -T.S/  is  n  xk.  The  number  of 
2  s  4  g  r 

samples  k  is  chosen  such  that  (-T^S/  -T.V  )  =  W  has  rank 

s  g 

nr  and  in  general, 

(157)  TjR'1  >  (T1Vs-CsYs  +  HcdUs)W-1  =  Z 

Again,  using  the  diagonal  structure  of  R,  (157)  can 
be  rewritten  as. 


(158)  (  r^t21,  ...  ,  rn  r— 2n  r  ^  =  1  '  ’•*  ' 


nr— 2nr  ‘ 


-nr' 


(158)  is  solved  for  the  r^,  i  =>  1,  2,  ...  ,  n^  by 
forming  the  column  sum. 
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(159) 


•  *  » 


n 


n 


£  fc2ki  ~  £  j  3  1#2,  . 

:=1  k=l 


'3k 


n 


(160) 


and  if  there  are  p  non-zero  column  sums  in  z, 
n  c 

k^i 11 2kj 

r_.  -  — — -  j  *  1,  2,  ...  ,  p  <  nr 


J ^ 


The  remaining  nr  -  p  diagonal  elements  of  r  must  be 
defined  by  additional  measurements.  A  similar  technique 
may  be  used  if  jL  is  taken  as  one  of  the  independent  vari¬ 
ables  instead  of  v  . 

~9 


ALTERNATIVE  SOLUTION  TECHNIQUE  FOR  R,  G 

An  alternative  technique  to  solve  for  the  diagonal 
elements  of  R,G  is  to  write  (120) ,  (121) ,  (122)  as  a  non¬ 
linear  set  of  algebraic  equations  in  the  form, 

2  2 

(161)  fi  *  *  "rnr#gl' *  '  ’  ,gng^"hi=0  i=1'2'  •  ♦  •  >n  +n i+ncn i 
and  solve  using  the  technique  of  successive  approximations 
[Rail,  M~7].  The  algorithm  is  an  iterative  one  given  by, 

(162)  xk+i  =  xk-J_1(xk) (f (xk)-h) 

and  x  is  the  vector  with  components  of  (161),  ff  (2£) ,  h  are 

)c 

the  vector  equivalents  of  (161)  and  J(x  )  is  the  Jacobian 
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matrix  of  (161)  evaluated  at  the  k-th  iteration.  The  tech¬ 
nique  of  successive  approximations  works  well  for  the  non¬ 
linear  equations  considered  because  the  types  of  non-line¬ 
arities  are  product  and  square  law  types. 

The  initial  starting  value  is  obtained  from  the  ap¬ 
proximate  solution  technique  presented  earlier.  The  major 
difficulty  is  the  algebra  involved  in  writing  out  the  equa¬ 
tions,  writing  out  the  form  of  the  Jacobian  matrix,  and 
eliminating  the  redundant  equations  so  that  J~^  exists. 

DETERMINATION  OF  R,G  -  SPECIAL  CASE  (COMPLETE  CIRCUITS) 

For  the  case  where  the  network-like  system  is  com¬ 
plete  in  the  sense  of  Bray ton  and  Moser  [FE-3J  all  of  the 

matrices  C  »  L  ,  R,  G  may  be  determined  without  additional 
s  s 

measurements . 

A  complete  circuit  is  a  network  which  has  special 
properties.  The  reader  is  referred  to  Brayton  and  Moser 
[FE-3,  p.  4  and  13]  for  a  thorough  discussion  of  complete 
circuits.  The  essence  of  the  completeness  property  is  that 
the  state  equations  can  be  defined  in  terms  of  a  potential 
function.  A  procedure  similar  to  that  for  formulating 
dynamical  equations  from  the  Hamiltonian  function  in  clas¬ 
sical  dynamics  can  be  employed  to  write  the  state  equations. 
For  a  complete  circuit,  the  state  equations  have  the 

form. 
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di 


dP 


(163)  Lj“dt  “  al.  j  1'2'  ***  '  nj e 

dv 

(164)  —§  k-  1.2 . nc 

JC 


For  the  linear,  constant  coefficient  case  P  has  the 


form. 


i  n 
1  r 


n 


(165)  P  =  2  zr-  ^  +  |  S\v%  -Ibjnvn+i^T  V. 

j=l  J  k=l  h=l  n=l 


where  i.,  j  =  1,2,  ...  ,  n  <n,v,  h  =  1,2,  ...  ,  n  <  n 
1  r  n  g  c 

are  components  of  the  state  vector the  e^'s  and  j^'s 

are  the  sources,  and  the  r.'s  and  g.'s  are  the  resistors 

3  ^h 

and  conductors. 

Applying  (163),  (164)  to  (165),  the  parameter  sub- 
matrices  have  the  forms. 


\ 


d66)  CsAn  = 


ng. 


*0 


(16?)  CsA12  -  Tx 
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0 


(168)  Lsa22  = 


'0 


(169)  LaA21  =  -T l 


where  A,  .  is  n  xn  and  A„„ 
li  c  c  22 

diagonal . 


is  n  xn ,  and  C  , 
l  l  s' 


L  are  both 
s 


(170)  CaBl  =  Hx 

(171)  LsB2  =  H2 


The  simplicity  of  the  above  formulation  allows  the 
determination  of  C  ,  L  ,  R,  G  by  inspection. 

5  S 


NUMERICAL  EXAMPLE 

Several  example  problems  are  presented  illustrating 
the  techniques  discussed.  Both  complete  and  other  types  of 
circuits  are  considered. 

Example  Problem  #1,  Second  Order  Bridge  Circuit 

A  bridge  circuit  has  some  features  which  make  it  a 
good  vehicle  for  illustration.  The  most  important  feature 
for  parameter  estimation  is  that  there  is  extensive  tree- 
cotree  resistive  coupling.  An  immediate  effect  of  the  tree- 
cotree  resistive  coupling  is  that  knowledge  of  che  system 
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states  and  topology  is  not  sufficient  to  determine  all 
parameters  uniquely. 

Figure  6  shows  the  circuit  schematic  and  its  associ¬ 
ated  graph.  The  state  variables  are  the  inductor  current, 
i^,  and  the  capacitor  voltage,  e^.  Equations  (172)  to 
(182)  define  the  circuit  model  and  parameter  matrices. 

The  main  solution  technique  using  linear  equations 
and  the  non-linear  technique  using  Newton- Raph son  iteration 
are  illustrated. 


(172)  A  =' 


/ 1  , 

1  1 

+  1 

1 

/  R1  ^ 

R3  1 

\ 

L  388 

.  056\ 

(csl 

iVR2 

R3+R4i 

C5 

1R1+R2  '  VR4 i 

u 

f 

i  I 

RX 

R3  ] 

I  | 

R1  r2  , 

R3  R4 ' 

n 

/ 

L  033 

-.570' 

\V 

^VR2 

R3+R4j 

L6 

1  Rl+R2  R3+R4- 

0  \ 


073) 

b  = 

1  i- 

Z2 

\  l6 

(174) 

T1  ° 

0 

(175) 

T2  = 

(  -1 

1) 

-1 

’ 

(176) 

T3  - 

;  -1 

1  0 

(177) 

t4  = 

,0  1 

/0  \ 

1-400/ 
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(178) 


R  = 


2.000  0  ' 
0  3.000 

J 


(179) 


(180) 

(181.) 

(182) 


K  °' 

i 

i 

a. ooo 

0  ' 

\  0 

R3 

=  | 

i  o 

1.000  , 

C  =  C,  =  1.500 
S  5 

L  =  Lr  =  2.500 


There  are  five  equations  and  six  unknowns,  therefore, 
one  of  the  parameters  must  be  specified.  If  a  unique  solu¬ 
tion  is  desired  then  additional  measurements  are  necessary 
for  unique  specification  of  parameters. 

Alternatively,  an  arbitrary  (positive)  value  may  be 
assigned.  If  the  latter  course  is  followed  then  the  others 
will  be  determined  by  that  specification. 

The  equations  to  be  solved  are  (183)  to  (187) . 


(183) 


b2  =  2.50 


(184) 


R, 


R. 


\ R2+Rl 


R3+R4 / 


“C5  *  a12 
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(185) 


R1R2 

Rl+R2 


R3R4 

R2+R4 


“  L6a22 


<186>  - 


R-+R . 
2  4 


-  c5alx 


(187) 


R, 


Rl+R2  R.3+R4 


L6a21 


Since  the  off-diagonal  terms  are  related,  it  follows 
that  may  be  computed  from  A,  b  alone. 

(188)  Cr  «  —  *  L,  =  1.50 

*  a12  8 


For  determining  the  resistive  values  there  are  three 
equations  and  four  unknowns,  thus,  a  resistive  value  must 
be  measured  or  specified. 

For  the  purpose  of  the  example,  it  was  assumed  that 
was  known  and  the  technique  of  successive  approximations 
used  to  solve  the  resulting  equations.  By  manipulation, 
the  non-linear  equations  to  be  solved  are. 


(189) 


R1R2 

VR2 


3R- 


.+3 


1.42 


(190) 


VR2 


R3+R4 


.58 


(191) 


.08 


Rl+R2 


V3 


To  solve  the  above  equations,  a  parameter  vector 


I  Ri 


r: 

below. 

(192) 


was  defined  and  used  in  the  Newton-Raphson  scheme 


K+l  k  -1 ,  kv  ,  k.  ,  . 
£  =  £  -  J  (£  )  U(£  )  -h) 


The  Jacobian  matrix  was, 


(1-R^)  (1-R*) 


\ 


12  13 


r5+3 


V3 


-  1 


3k 

V 


3 


Ikil)  j 

(R3+3)/ 


An  initial  guess  which  was  50  percent  in  error  was 
used  to  start  the  Newton-Raphson  iteration.  The  results 
are  shown  in  Table  7. 
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Table  7.  Newton- Raphson  Iterative  Solution  for 
Parameter  Values 


PARAMETERS 

ITERATION  NUMBER 

Rl 

R2 

R3 

0 

1.5 

1.0 

1.5 

1 

1.62 

1.35 

0.76 

2 

0.96 

2,23 

0.68 

3 

0.61 

2.16 

1.20 

4 

* 

0.82 

2.11 

1.13 

5 

0.96 

2.03 

1.01 

6 

0.99 

2.01 

1.01 

7 

1.00 

2.00 

1.00 

Correct  Values 

1.00 

2.00 

1.00 

To  illustrate  the 

main  solution 

technique 

it  was  as- 

sumed  that  an  identification  scheme 

had  yielded 

the  A,  b 

matrices  and  the  state 

trajectories. 

The  state 

derivative 

was  computed  by  x.  s  Ax  +  b  u.  An  additional  voltage  measure¬ 


ment  was  required  because  the  rank  of  the  matrix  T-j  was  1. 

To  determine  the  vector  v  =  (  1)  uniquely,  the  relation 

c3 

(150)  is. 

(194)  -T,  **v  «*  Li  +  u 

3  — g  s— s 


64 


Substituting  for  T^,  i^,  Lg,  gives  the  branch  volt¬ 

age  relation  below. 

(195)  -(e1  +  e3)  =  -L&i6  +  e 

By  measurement  of  e^,  e^  can  be  computed  since  L^, 
e,  are  known. 

Equation  (156)  was  used  to  determine  the  chord  resis¬ 
tors  r2<  r^. 


(196) 


(_A_  -L, 

r2  V 


e5(ti)-e1(t1)e5(t2)-e1(t2) 
— e^  (t^)  — e^  — ^5  ^^”1^  “e3 


c5{e5(tl)^5{t2)' 


Samples  were  taken  at  a  1.0  sec,  t2  =>  3.2  sec. 
Table  8  shows  the  sampled  values  and  the  results  for  r2# 

V 
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SECOND  ORDER  BRIDGE  CIRCUIT  STATE  AND  STATE  DERIVATIVE 
RESPONSE  TO  STEP  OF  TEN  VOLTS 


MEASURED  AND  COMPUTED  BRANCH  VOLTAGE  FOR 
SECOND  ORDER  BRIDGE  CIRCUIT 


O 


so&e^ioA  tputua 


oo 


8  V 


FIGURE 


Table  8.  Computational 

Results  for  Chord 

Resistors 

VARIABLE 

TIME 

t  =  1.0  sec. 

t  =  3.2  sec. 

cb®5  “  ^P® 

-.208 

-.235 

e^  -  Volts 

-.082 

-.444 

e  -  Volts 

10 

10 

e^  -  Volts 

-2.059 

-4.063 

e^  -  Volts 

-2.266 

-4.293 

r2  -  2  ohms 

- 

- 

-  3,01  ohms 

Correct  Values 

r2  53  2  ohms 

52  3  ohms 
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Once  the  chord  resistors  were  known, the  branch  resis¬ 
tors  were  easily  determined  by. 


(197) 

rl  =  el  (fcl)  ^1  ) 

(198) 

r3  =  e3(tl)(i3(tl}) 

The  branch  currents  were  given  by  the 

chord  currents 

where , 

(199) 

h  =  'V  i2  ■  _i6  +  e2/v2  ”  i6  + 

(v^) 

r2 

(200) 

h  "  _i6  +  H  H  *  Vv4  "  -'e  + 

(— ec— e^) 

r4 

At 

t^  =  1.0  sec,  the  values  of  i^,  i^ 

were  computed 

below. 

(201) 

(1.0)  «  1.0 

(202) 

i 3 (1.0)  =■  1.0 

The  result  of  the  example  is  that  given  two  state 
measurements,  an  additional  branch  voltage  measurement, 
and  the  input,  that  all  six  of  the  network  parameters  can 
be  computed. 
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Example  Problem  #2  -  Sixth  Order,  Two-input  Complete 
Circuit. 

Example  problem  #2  illustrates  the  parameter  determina¬ 
tion  technique  for  the  special  case  of  a  complete  circuit. 

The  system  is  the  complete  circuit  for  Example  #2.  Table 
9  shows  the  results  of  the  numerical  computation  and  the 
corresponding  parameters. 

Some  of  the  parameters  are  over-determined,  that  is, 
there  are  more  than  one  equation  for  some  parameters.  A 
useful  benefit  of  the  above  fact  is  that  computational 
"noise"  can  be  smoo tiled  out  by  averaging  results.  A  check 
on  the  consistency  of  the  identification  is  also  provided. 

Table  10  shows  the  results  of  the  computation  for 
each  parameter  when  the  average  is  used  for  the  over-de¬ 
termined  equations. 

Example  problem  2  shows  that  measurement  of  six  states 
and  two  inputs  allows  the  computation  of  II  parameters. 

The  computation  involves  simple  algebraic  techniques.  Ex¬ 
ample  problem  #1  contrasts  with  problem  #2  in  the  relatively 
complex  technique  of  computation  plus  the  need  for  additional 
measurements  other  than  state  and  input. 


Table  9.  Numerical  Results  for  Complete  Circuit. 


PARAMETERS 

COMPUTED  VALUES  FROM 
IDENTIFICATION 

i/r4Cj 

5.01 

1/C1 

10.00 

1/R5c2 

2.00 

1/C2 

2.00 

1/C3 

1.00 

l/c3 

1.00 

1/C3 

1.00 

i/l*! 

10.02 

1/1L1 

9.96 

l/L1 

10.02 

VLi 

4.00 

i/l2 

5.00 

i/l2 

5  >'r 

1A2 

5.00 

VL2 

5.00 

1/Lj 

3.34 

IA3 

3.34 

R3/L3 

10.00 
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Table  10.  Comparison  of  Actual  and  Computed  Parameters 
For  Complete  Circuit. 


PARAMETER  COMPUTED  VALUE  ACTUAL  VALUE 


C1 

0.10 

0.10 

C2 

0.50 

0.50 

C3 

1.00 

1.00 

R1 

1.01 

1.00 

r2 

1.00 

1.00 

R3 

3.00 

3.00 

R4 

1.99 

2.00 

R5 

1.00 

1.00 

L1 

0.10 

0.10 

L2 

0.20 

0.20 

L3 

0.30 

0.30 

Example  Problem  #3  -  Complete  Circuit  Model  of  Second- 
Order  Bridge  Circuit. 

A  pertinent  issue  for  modeling  purposes  is  the  exist¬ 
ence  of  a  model  of  the  system  which  will  give  the  same 
response  as  that  measured  but  have  a  different  topology 
and  parameter  values.  Example  #3  shows  that  a  model  of  the 
circuit  in  Example  #1  exists, which  has  a  different  topology 
and  different  parameter  values,  yet  has  the  same  A-matrix. 
The  input  matrices  are  different,  thus,  for  the  same  input, 
different  state  responses  will  be  observed. 

A  multiplicative  scaling  of  the  bridge  circuit  input, 
however,  to  provide  the  input  to  a  complete  circuit  model 
of  the  incomplete  circuit  yields  the  identical  state  tra¬ 
jectory. 

Tat  net  rtsult  is  that  there  exists  a  second-order 
complete  circuit  model  of  the  second-order  bridge  circuit 
which  has  identical  stability  properties,  i.e.,  has  the 
same  A-matrix  and  the  two  circuits  have  identical  free 
responses  to  identical  initial  conditions. 

Additionally,  the  forced  responses  to  the  same  input 
differ  only  by  a  multiplicative  constant.  Figure  9  shows 
a  schematic  of  each  circuit. 
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FIGURE 


The  A,  b  matrices  of  the  bridge  circuit  are. 


(203)  A  = 


-.388 

.033 


-.056 

-.570/ 


(Bridge  Circuit) 


(204)  b  = 


.4 


(Bridge  Circuit) 


By  comparison,  the  A,  b  of  the  complete  circuit  are. 


-1 


-1 


(205)  A  = 


R1C2 


2 

-R 


4 


(Complete  Circuit) 


(206)  b 


0  \ 


\L3  / 


(Complete  Circuit) 


If  the  A,  b  matrices  were  identical,  then  an  incon¬ 
sistency  would  arise  since  there  are  inconsistent  require¬ 
ments  for  L-j,  i .  e . ,  from  the  a21'  ^2  elements. 

(207)  L,  -  — ~  »  30  h 

a21 

(208)  L3  -  l/b5  »  2.5  h 

By  choosing  the  complete  circuit  parameters  to  satisfy 
(consistently) , 
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!W*: 


Oley&WW 


(209)  C2  =  17.90  farads 

(210)  =  6.95  ohms 

(211)  =>  30  henries 

(212)  R5  =  (30)  (.570)  ohms 

then,  the  complete  circuit  has  the  same  A-matrix  as  the 
bridge  circuit.  Additionally,  if  the  input  e(t)  is  multi¬ 
plied  by  a  constant,  k,  and  used  to  drive  the  complete  cir 
cuit,  then  identical  state  trajectories  for  arbitrary 
initial  conditions  result.  The  gain  constant  must  satisfy 

(213)  k  •  .4L3 

Alternatively,  it  may  be  desired  to  choose  (the 
only  inconsistently  specified  parameter)  to  minimize  some 
error  criterion. 

Referring  to  the  results  of  Chapter  ~  the  state 
error  equation  is. 


!  -  .380 

-.0561 

f  0 

0  \ 

(214:  £  "  S  1/Lj 

-.570 ) 

e  +•  | 

~  ! 

■ 

.033-1 

1  I-3 

1 

0 

/ 

0  \  A 

-  +  4-1  i  U'“  S  - 
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The  A,  b  matrices  of  the  bridge  circuit  are. 


(203)  A  * 


I-.388 


.033 


-.056 

-.570, 


(204)  b  * 


0  ’ 

.4 


(Bridge  Circuit) 


(Bridge  Circuit) 


By  comparison,  the  A,  b  of  the  complete  circuit  are, 


/  r'c 


(205)  A  = 


12 


2 

-R , 


(Complete  Circuit) 


(206)  b 


-/ 


0  \ 


\L3 


(Complete  Circuit) 


If  the  A,  b  matrices  were  identical,  then  an  incon¬ 
sistency  would  arise  since  there  aro  inconsistent  require¬ 
ments  for  L^,  i ♦ e . ,  from  the  a.,^,  ^2  elements. 


(207)  L  j  »  —  =  30  h 

J  a21 

(208)  L3  =  l/b2  «  2;  5  h 

By  choosing  the  complete  circuit  parameters  to  satisfy 
(consistently) , 
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5fH  * 


(209)  C2 

(210)  ^ 

(211)  L3 

(212)  Rg 


17  ..90  farads 
6.95  ohms 
30  henries 
(30)  (.570)  ohms 


then,  the  complete  circuit  has  the  same  A-raatrix  as  the 
bridge  circuit.  Additionally,  if  the  input  e(t)  is  multi¬ 
plied  by  a  constant,  k,  and  used  to  drive  the  complete  cir¬ 
cuit,  then  identical  state  trajectories  for  arbitrary 
initial  conditions  result.  The  gain  constant  must  satisfy, 

(213)  k  =  .4L, 


Alternatively,  it  may  be  desired  to  choose  L3  (the 


only  inconsistently  specified  parameter)  to  minimize  some 
error  criterion. 


Referring  to  the  results  of  Chapter  3,  the  state 
error  equation  is. 


i 

-.388 

-.056 

!  0 

0  \ 

(2b*)  e  » 

i 

i/l3 

-.570 

£  + 

.033-1 

l  l3 

0 

10  \  A 

X  +| .4-1  U'“  =  - 

l7/ 
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One  criterion  which  could  be  used  is  to  minimize  the 
mean  square  steady  state  error  to  a  step  and  require  L3>  0. 
Solving  for  the  steady  state  error  (u  -  10  volts,  x^  =  -1.0 
volts  at  steady  state) . 


(215)  el/ss  -  -144e2ss 


(216)  e2/ss  “  (1.024) (L3  3.967-9) 

IT  5  7  0L  3  - .  14  ■* 

The  mean  sq*  ->rt  error  (  e*}  is  the  sum  of  the  squares 
of  the  mean  error  components. 


(217)  e  e1/gs  +  e3/gs  (1*024)i-1.570-.144 


(L33. 967-9.0) 


A  minimum,  e  =*  0  exists  for. 


(218)  L3  (3.967)  =  9  or 
L3  =  2.27 

Figure  10  shows  the  corresponding  error  trajectory. 
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FIGURE  10 


SUMMARY 


Chapter  7  has  dealt  with  techniques  for  extracting 
the  parameter  submatrices  from  the  A,  B,  matrices.  System 
tooolony  plays  an  important  role  because  its  knowledge 
enables  one  to  determine  a  priori  if  knowledge  of  the  state 
and  input  is  sufficient  to  determine  all  parameters  uniquely. 
Such  knowledge  is  mportant  in  planning  a  testing  program. 

If  one  tries  to  determine  the  parameters  from  know¬ 
ledge  of  the  A,  B,  matrices  anc'  system  structure  on./,  then 
a  set  of  nonlinear  algebraic  equations  must  be  solved. 

Newton -Raphson  iteration  is  eff active  for  this  case. 

For  the  case  where  additional  tree  or  cotree  variables 
can  be  measured,  then  the  parameters  c*n  be  determined  by 
linear  algebraic  operations. 

Example  problems  illustrating  the  various  techniques 
were  presented. 
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APPENDIX 

FORMULATION  OF  NE WORK-EQUATIONS  IN  STATE  FORM 

The  following  development  is  based  on  the  work  by 
MacFarlane  [F2-9]  with  modifications  in  notation  and  form 
to  suit  the  problem  of  parameter  estimation. 

Let  Tq  be  the  matrix  relating  branch  and  chord  vari¬ 
ables  for  a  tree  of  the  network  which  contains  the  maximal 
number  of  capacitive  elements  and  such  that  the  cotree  con¬ 
tains  the  maximal  number  of  inductive  elements.  The  only 
capacitive  elements  not  in  the  tree  are  those  which  form 
capacitor  only  loops  with  capacitors  in  the  tree,  and  the 
only  inductive  elements  which  are  not  in  the  cotree  are 
those  which  form  inductor  only  cut  sets  with  inductors  in 
the  cotree.  Also,  the  tree  is  selected  to  include  all 
voltage  sources  such  that  all  current  sources  are  in  the 
cotree.  If  the  currents  ana  voltages  are  selected  as  above 
then  the  vector-matrix  equations  given  below  follow. 


'  ^  1 

1  0  Tol 

1% 

(219  ) 

t  ~c  ) 

sa 

K  °  i 

U  J 

where  the  subscripts  b,  c  denote  branch  and  chord  variables 
respectively. 
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The  submatrix  given  by  (223)  is  the  matrix  for  the 
case  where  all  voltage  sources  are  open-circuited  and  all 
current  sources  are  short  circuited. 


(223) 


The  component  equations  are  assumed  to  be  of  the  form# 


(224a) 

(224b) 

(224c) 

(224d) 


ibc  =  Gs4jc 
^bg  =  G  ^bg 


/rsi' 

f^s_ 

<} 

1 

iMi 

Le  J 

i  ^b/  i 

The  state  variables  are  defined  to  be# 


A 

(225a) 

(225b)  4^, 


i 

i 

I 

i 


[ 

i 


) 
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!’■ }  ^  ffW.1  ■  p  V  5  %-| 


T*'*’  ^7,L,irp-A  ^i;^- w.yr'*r*-«~TT'  ••.  -*vr>  '•,  --1--^"  .v.'.rav^  ^vi 


The  remaining  variables  are  renamed# 


(226a) 

±br  -  it 

(226b) 

hat  a  ie 

(226c) 

^bg  =  ¥g 

(226d) 

£cc  =  2e 

The  subscripts  s,  r,  g,  e,  denote  state,  resistive, 
conductive  and  excess  dynamic  element  vector,  respectively. 

By  substituting  the  renamed  variables  and  component 
equations  into  (219)  and  rearranging,  the  network  equations 
are  in  (225) . 


Cs^s 


l  1+mH 

s— s  x— e 


(227) 


Gv 

“9 


^r 


0  Tl  0  t2 


i 


-T^O  -T3O 


0  T3  0  t4 


-t2o  -T40 


-t£o  0  0 


,0  t5  0  0 


Leie+Mi U 


H 


‘cd 


H 


gd 


id 


^cd 


^bd 


rd 


K 


cd/ 


'll 
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f 

V 

y 

where, 

C  diagonal  n  xn  matrix  of  tree  capacitors 

S  CO 

Cg  -»  diagonal  n^xn^  matrix  of  chord  capacitors 

L  -♦  symmetric,  positive  definite  n  xn  matrix 
s  t  i 

of  chord  inductors. 

Le  -»  symmetric,  positive  definite,  n^n^  matrix 
of  tree  inductors. 

G  diagonal  n  xn  matrix  of  tree  conductors 
9  9 

R  diagonal  nrxn  matrix  of  chord  resistors 

From  the  dimension  definitions  of  the  component  raatri 
ces,  the  topological  matrices  have  dimensions. 


T.  -♦  n  xn  . 
1  c  X 


T0  -♦  n  xn„ 
2  c  r 


T3  -  ngxnl 
T4  -  ngxnr 


Tc  n„xn. 
b  2  l 


TV  -♦  n  xn, 
o  cl 


The  input  matrices  have  dimensions, 

Kjd  -  vpc 

Krd  -  nrxpc 
Kcd  -  nlxPc 
Hcd  -  noxpb 
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V  " 

Hid  -  n2XPb 


By  eliminating  the  excess  variables ,  the  equation  has 
the  form, 


(228)  Jx  =  Fx  +  Hu 
where, 


(229) 


(230) 


(231) 


(232) 


(233) 


C  +TcC  TV 
s  6  e  6 


VT^i  +  T^LeT5J 


H 

"Vh 

E* 

o 

1 _ 

0  T2' 

|G  T4] 

1°  T3' 

•* 

Ll-£  •/ 

*1 

\-Tt  *1 

i-T£  °l 

H  = 


/'  0 


t2.  !g 


TV 


\"1 


-T3  °i  i  -T4  R 


'agd  0  \:'Hoa 


0  \ 


Krd/:l° 


K 


Id  I 


o 

o 

(T  6 

0  ' 

:/C 

•  1  w 

0  1 

1° 

Kcd' 

\-Mi  Hcd°l 

0 

-T^ 

A5/ 

*1  0 

LeJ 

lHid 

°l 

v  \ 

— s 


.  i 
\  — s 


^cd 


u  =  i 


^bd 


-cd 

^bd  I 
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In  summary  of  subscript  notation,  s  -*  state,  g  -♦  con¬ 
ductive,  r  -*  resistive,  e  -♦  excess  dynamic,  cd  -*  chord 
driver  (current  source) ,  and  bd  -»  branch  driver  (voltage 
source) . 
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